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Abstract The difference between the electromagnetic 
self-energies of proton and neutron can be calculated 
with the Cottingham formula, which expresses the self¬ 
energies as an integral over the electroproduction cross 
sections - provided the nucleon matrix elements of the 
current commutator do not contain a fixed pole. We 
show that, under the same proviso, the subtraction 
function occurring in the dispersive representation of 
the virtual Compton forward scattering amplitude is 
determined by the cross sections. The representation in 
particular leads to a parameter-free sum rule for the nu¬ 
cleon polarizabilities. We evaluate the sum rule for the 
difference between the electric polarizabilities of proton 
and neutron by means of the available parameteriza- 
tions of the data and compare the result with experi¬ 
ment. 
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1 Introduction 

The mass difference between proton and neutron had 
been puzzling for a long time. Ever since Heisenberg had 
introduced isospin symmetry to explain the near degen¬ 
eracy of these two levels [T], it was taken for granted 
that the strong interaction is invariant under isospin 
rotations and that the mass difference is of electromag¬ 
netic origin. In this framework, it was difficult, however, 
to understand the experimental fact that the neutral 
particle is heavier than the charged one. A first step 
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towards a resolution of the paradox was taken by Cole¬ 
man and Glashow, who introduced the tadpole domi¬ 
nance hypothesis Hi, which associates the bulk of the 
electromagnetic self-energies with an octet operator. 
The origin of the tadpole remained mysterious, how¬ 
ever. The puzzle was solved only in 1975, when it was 
realized that the strong interaction does not conserve 
isospin, because the masses of the up- and down-quarks 
strongly differ [4] . The crude estimates for the ratios of 
the three lightest quark masses obtained in that work, 
m„/md — 0.67, m s /md — 22.5, have in the meantime 
been improved considerably. In particular, Weinberg j5j 
pointed out that in the chiral limit, the Dashen theorem 
provides an independent estimate of the quark mass ra¬ 
tios, as it determines the electromagnetic self-energies 
of the kaons in terms of those of the pions. Neglecting 
higher orders in the expansion in powers of m u ,md , 
and m s , he obtained the estimate m u /rrid — 0.56, 
m s /md — 20.1. Also, the decay 77 —>- 37 t turned out 
to be a very sensitive probe of isospin breaking [SHU!- 
The quark mass ratios obtained from that source also 
confirmed the picture. According to the most recent 
edition of the FLAG review E3, the current lattice av¬ 
erages are m u /md = 0.46(3), m s /md = 20.0(5). 


1.1 Cottingham formula, dispersion relations 


While the functions V\{v, q 2 ), V 2 (is, q 2 ) are odd 
under v —> —v, the time-ordered amplitudes 
Ti{v,q 2 ),T 2 {y,q 2 ) are even. In view of the contribu¬ 
tions arising from Regge exchange, V\(y,q 2 ) ~ v a , 
V 2 (v,q 2 ) ~ only T 2 obeys an unsubtracted dis¬ 

persion relation, while for Tj a subtraction is neededQ 
For q 2 < 0, the dispersion relations thus take the form 


T 1 (u,q 2 )=S 1 (q 2 ) + 2u 2 

pOO 

T 2 (u, q 2 ) = 2 / dv' v' - 

Jo 1 


[°° dv' Vi(y', q 2 ) 

/q // v' 2 — v 2 — ie ’ 

V2W,q 2 ) 

— 7/2 — 7 V 


( 2 ) 


The formulae hold in the cut ^-plane; the upper and 
lower half-planes are glued together along the inter¬ 
val \v\ < Q 2 /2m of the real axis (throughout, we 
use Q 2 = —q 2 whenever this is convenient). As illus¬ 
trated with the discussion in appendix 11 it is impor¬ 
tant that kinematic singularities, zeros, and constraints 
be avoided - throughout this paper, we work with the 
amplitudes defined in appendix [A] which are free of 
these mm- 

We refer to S±(q 2 ) as the subtraction function. It 
represents the value of the amplitude T\(y , q 2 ) at v = 0. 
For later use we introduce the analogous notation also 
for T 2 (y , q 2 ): 


Si{q 2 ) = Ti(0, q 2 ), S 2 (q 2 ) = T 2 (0, q 2 ). (3) 


The analysis of [4] relies on the Cottingham for¬ 
mula [ T2] . which invokes dispersion relations to re¬ 
late the spin-averaged nucleon matrix elements of 
the time-ordered product, {p\Tj ,J ’(x)j v (y)\p), to those 
of the commutator of the electromagnetic current, 
(p\[j fJ '{x),j ,y (y)]\p). Lorentz invariance and current con¬ 
servation determine the Fourier transforms of these 
matrix elements in terms of two invariant amplitudes, 
which only depend on the two variables v = p ■ q/m 
and q 2 , where m is the nucleon mass and q the photon 
momentum. We stick to the notation used in [4] and 
denote the invariant amplitudes by T\{y,q 2 ),T 2 {y : q 2 ) 
and Vf (v, q 2 ), V 2 (v, q 2 ), respectively. Explicit formulae 
that specify the matrix elements (p\Tj li (x)j 1 '(y)\p) and 
j"(y))\p) i n terms of the invariant ampli¬ 
tudes are listed in appendix E where we also exhibit 
the relations between the structure functions V\(v, q 2 ), 
V 2 (v,q 2 ) and the cross sections ctt and ox of electron 
scattering. 

In the space-like region and for v > 0, the structure 
functions represent the imaginary parts of the time- 
ordered amplitudes: 

Im T\(y, q 2 ) = nVi(v,q 2 ) , 

ImT 2 (i/, q 2 ) = irV 2 (v,q 2 ) , v > 0 , q 2 < 0 . (1) 


1.2 Reggeons and fixed poles 


In [1] it is assumed that the asymptotic behaviour is 
determined by Reggeon exchange. The contribution of 
a Regge pole to a scattering amplitude at large center- 
of-mass energy squared s and small momentum transfer 
t < 0 has the form (see e.g. my- 

T(s,t) = - nl3 ° ^ {exp[—?7ra(t)] + r}s a(t) , (4) 

sin7ra(t) 


where a(t ) and /3(t) denote the trajectory and the 
residue, respectively, and r is the signature. In the con¬ 
text of the present paper, we are concerned with t = 0 
and t = 1. The continuation of the asymptotic for¬ 
mula (J4]) to low energies is not unique. For definiteness, 
we work with the representation 




sin ira 

a>0 


x | (so - s+ - ie) a + (s 0 - s_ - ie) a j , (5) 


1 It was even suggested that the subtraction term might solve 
the notorious puzzle with the proton-neutron mass difference, 
as it could dominate over the remainder and explain why the 
neutron is heavier than the proton |13) . In hindsight, we know 
better. 
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where s+ and s_ stand for s± = (p ± q ) 2 3 = m 2 ± 
2mi'—Q 2 and So > m 2 is a constant. The expression m 
is manifestly symmetric under photon crossing. Unless 
the intercept a is an integer^ the first term in the curly 
brackets contains a branch cut along the positive real 
axis, starting at 2 mv = So~m 2 + Q 2 . The second is real 
there. One readily checks that, on the upper rim of this 
cut, the individual terms in the sum 0 differ from the 
asymptotic expression 0 only through contributions 
of 0(s“ -1 ). 

The basic assumption made in 2] is that, in the 
limit v —> oo at fixed q 2 , only the Reggeons survive, so 
that the difference tends to zeroU 

Ti(i/,g 2 )-?«(!/, 9 2 )->- 0, (6) 

We refer to this hypothesis as Reggeon dominance. 
A nonzero limiting value in (0 would represent a v- 
independent term. In Regge-language, a term of this 
type would correspond to a fixed pole at angular mo¬ 
mentum J = 0. The Reggeon dominance hypothesis 0 
thus excludes the occurrence of such a fixed pole. 

The presence or absence of a fixed pole at J = 0 in 
Compton scattering is a standard topic in Regge pole 
theory [T8] and the literature contains several works 
advocating the presence of such a contribution. In par¬ 
ticular, the universality conjecture formulated in m 
has received considerable attention (see e.g. m and 
the papers quoted therein). 

Note, however, that these considerations go be¬ 
yond the safe grounds provided by asymptotic freedom. 
While the short distance properties of QCD ensure that, 
if both v and q 2 are large, the behaviour of Tj (v, q 2 ) and 
T2{y,q 2 ) is governed by the perturbative expansion in 
powers of the strong coupling constant, the behaviour 
in the Regge region, where only v becomes large, is not 
controlled by the short distance properties of QCD. In 
particular, values of q 2 of the order of Aq CD are outside 
the reach of perturbation theory, even if v is large. 

The perturbative analysis shows that an infinite set 
of graphs needs to be summed up to understand the 
high-energy behaviour of the amplitudes in the Regge 
region. The dominating contributions can be repre¬ 
sented in terms of poles and cuts in the angular mo¬ 
mentum plane (Reggeon calculus, Reggeon field the¬ 
ory). The behaviour of the sum thus differs qualitatively 
from the one of the individual diagrams. 

There is solid experimental evidence for the presence 
of Reggeons also in the data. The relation 0 amounts 

2 Integer values of a require special treatment, but since this 
case does not arise for the parameterizations we are working 
with, we do not discuss it further. 

3 More precisely, it is assumed that the difference disappears 
rapidly, so that it obeys an unsubtracted dispersion relation. 


to the assumption that the asymptotic behaviour of the 
current correlation function can be understood in terms 
of these. In the analysis described in the present paper, 
this assumption plays a key role. In particular, as will 
be demonstrated explicitly below, it uniquely fixes the 
subtraction function relevant for the difference between 
proton and neutron in terms of the electron cross sec¬ 
tions, so that the entire self-energy difference can be ex¬ 
pressed in terms of these cross sections. In other words, 
the necessity of a subtraction in the fixed-*/ 2 dispersion 
relation for Ti(v,q 2 ) modifies the relation between the 
self-energy difference and the electron cross sections, 
but does not destroy it0 

The subtraction functions occurring in the fixed-t 
dispersion relations relevant for real Compton scatter¬ 
ing are analyzed in mm- As shown there, the experi¬ 
mental information on the differential cross sections can 
be used to impose bounds on the subtraction functions. 
In particular, these bounds lead to the conclusion that 
the electric polarizability of the proton is necessarily 
larger than the magnetic one, in conformity with exper¬ 
iment. An update of this work with the data available 
today is highly desirable. Unfortunately, this approach 
to the problem cannot readily be extended to virtual 
Compton scattering, because data on the differential 
cross sections are available only for real photons. 


1.3 Recent work 

The numerical analysis of j3] was based on the scaling 
laws proposed by Bjorken [23j. The data available at the 
time were perfectly consistent with these, but Bjorken 
scaling correctly accounts for the short-distance proper¬ 
ties of QCD only to leading order in the perturbative ex¬ 
pansion in powers of a s . The higher-order contributions 
generate specific violations of Bjorken scaling [MUSS]. In 
the meantime, the implications of the phenomenon and 
the corresponding modification of the short-distance 
properties of the matrix elements {p\Tj IJ, (x)j L '(y)\p) 
have been investigated by Collins [26]. Unfortunately, 
however, he did not reevaluate the self-energy differ¬ 
ence in this framework. In fact, the question of whether 
the Reggeons do dominate the asymptotic behaviour or 
whether the amplitude in addition contains a fixed pole 
at I = 1, J = 0 is not touched at all in that work. 

Motivated in part by the study of hadron electro¬ 
magnetic mass shifts on the lattice (see, e.g., [2IH2S]), 

4 In U, this conclusion was derived on the basis of a somewhat 
weaker form of Reggeon dominance, which does not invoke 
the matrix elements of the time-ordered product, but those of 
the current commutator. For a brief discussion of this aspect, 
we refer to appendix Id 
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the Cottingham formula has recently been reexam¬ 
ined [HIMTB] . but the central issue in this context - the 
possible occurrence of fixed poles - is not addressed in 
these papers, either. Instead, the electron cross sections 
ctt, and the subtraction function S±(q 2 ) are treated 
as physically independent quantities. The main prob¬ 
lem with the framework set up in ED is that a direct 
experimental determination for Si(q 2 ) is not available. 
To bridge the gap, the authors set up a model which pa¬ 
rameterizes the dependence of the subtraction function 
on q 2 . The overall normalization, S'i(O), can in princi¬ 
ple be determined from the difference between the mag¬ 
netic polarizabilities of proton and neutron, albeit the 
experimental value is subject to rather large uncertain¬ 
ties [57;. The main problem in this approach, however, 
is the momentum-dependence of the subtraction func¬ 
tion, which leads to a systematic uncertainty that is 
difficult to quantify. 

1.4 Structure of the present paper 

The remaining sections are organized as follows. In sec¬ 
tion H we show how the Reggeon dominance hypothe¬ 
sis © fixes the subtraction function Si(q 2 ) from space- 
like data alone. In section [51 we discuss the splitting 
of the amplitudes Ti into elastic and inelastic contri¬ 
butions. We derive sum rules for the nucleon polariz¬ 
abilities in section [D while a thorough phenomenologi¬ 
cal analysis is provided in section [5] In particular, the 
sum rules allow us to predict the difference between 
the electric polarizabilities of proton and neutron. In 
view of the fact that the proton polarizabilities are ex¬ 
perimentally known more accurately, our result can be 
turned into a prediction of the electric polarizability 
of the neutron, which is consistent with observation 
but somewhat more precise. The magnetic polarizabil¬ 
ities then follow from the Baldin sum rule. Section ED 
is devoted to the electromagnetic self-energies of pro¬ 
ton and neutron. We discuss the renormalization of the 
Cottingham formula, in particular the role of the sub¬ 
traction function in the evaluation of the self-energy 
and provide a comparison with recent work on the is¬ 
sue. A summary and concluding remarks are given in 
section [7] In appendix lAl we detail the notation used. 
Appendix [Bl reviews those properties of Compton scat¬ 
tering we are making use of. In particular, we discuss 
the frame-dependence of the spin average and derive 
the low-energy theorem which underlies the sum rule 
for the electric polarizability. Appendices [C] and [D] con¬ 
tain a short discussion of the role of causality in our 
analysis. Last but not least, we note that in [5UH55j . a 
comparison with the analysis of [1] is attempted. Unfor¬ 
tunately, many of the statements made there are sim¬ 


ply incorrect. Some of the misconceptions are rectified 
in appendix [E] 


2 Determination of the subtraction function 


The Regge amplitude obeys a once-subtracted fixed-q 2 
dispersion relation: 


nOO 7 / 

T R (v,q 2 )=T R (0, q 2 ) + 2v 2 4 

Jo v 


dv' v/V,? 2 ) 


( 7 ) 


In the space-like region and for v > 0, the absorptive 
part of the amplitude specified in © is given by 


V R (v, q 2 ) = ^2 /3 a (Q 2 ) 0(s+ - s 0 ) (s+ - *>)“ • (8) 

a>0 

The Reggeon dominance hypothesis © implies that 
the difference between the full amplitude and the Regge 
contributions, T\{y,q 2 ) = Ti(V, q 2 ) — T R (v,q 2 ), obeys 
an unsubtracted dispersion relation. In particular, the 
value of Ti(0, q 2 ) = S\(q 2 ) — T R (0,q 2 ) is given by an 
integral over the difference q 2 ) — V R {v, q 2 ). Hence 
the subtraction function can be represented as 

S 1 (q 2 ) = T R (0,q 2 ) 

+ 2[°°-{V 1 (v,q 2 )-V R (v,q 2 )} . (9) 

Jo v 

This formula explicitly represents the subtraction func¬ 
tion in terms of measurable quantities: the structure 
function V\(v, q 2 ) is determined by the cross sections for 
inclusive electron-nucleon scattering. The high-energy 
behaviour of these cross sections also determines the 
Reggeon residues /3 a (Q 2 ) and thereby fixes the term 
T R (0,q 2 ), as well as the corresponding contribution to 
the structure function, V R {u,q 2 ). 

If the trajectory intercepts a were all below zero, the 
unsubtracted dispersion integral over V R (y, q 2 ) would 
converge and would exactly compensate the first term 
on the right of © the subtraction function would then 
be given by the unsubtracted dispersion integral over 
Vi (V, q 2 ). The expression for the subtraction function 
in © shows how the divergence of the unsubtracted 
dispersion integral generated by the Reggeons is han¬ 
dled: the corresponding contribution is removed from 
the integrand, so that the integral converges also at 
the physical values of the intercepts. The modification 
is compensated by the term T R (0,q 2 ), which must be 
added to the integral over the remainder. The procedure 
amounts to analytic continuation in a from negative 
values, where T R ( 0, q 2 ) is given by the unsubtracted 
dispersion integral over V R (u, q 2 ) to the physical val¬ 
ues, where that representation does not hold any more. 




5 


We emphasize that the specific form used for the 
Regge parameterization does not matter. In particular, 
the Regge amplitude specified in (0 involves a free pa¬ 
rameter, so- Since it does not affect the leading term 
in the asymptotic behaviour, the value used for so is 
irrelevant - our results are independent thereof. In the 
following, we simplify the equations by taking so in the 
range so > (m + M,) 2 , which has the advantage that 
V^iy.q 2 ) then vanishes outside the inelastic region. 


3 Elastic and inelastic contributions 

3.1 Elastic part 


In the case of T\{y, q 2 ) there are two contributions, one 
from the subtraction function, the other from the sub¬ 
tracted dispersion integral: 


/■= 

T?{v iq 2 ) = Sf{q 2 ) + 2 V 2 / 

Jo 


dv' V?(y',q 2 ) 


( 12 ) 


The sum rule 0 for the subtraction function implies 


Sf{q 2 ) = 2 T-Vf{u,q 2 ). (13) 

Jo v 

Taken together, the two terms on the right hand side 
of (fl2l) yield the unsubtracted dispersion integral, so 
that the expression takes the same form as the one for 
TZ\v,q 2 y. 


The contributions to the structure functions arising 
from the elastic reaction e + N —> e + N are determined 
by the electromagnetic form factors of the nucleon. In 
the space-likc region, these contributions are restricted 
to the lines q 2 = ±2 vm and read (i = 1, 2) 

Vi\v,q 2 ) = vf(q 2 ) {S(q 2 + 2 \mv) - S(q 2 - 2 mv)} , 

V l(l 2 ) = 4 m 2 _ q 2 { G b('? 2 ) - G m(<? 2 )} i 

el, 2 n _ 2m 2 ' 

V2[q (-q 2 )(4m 2 ~ q 2 )’ 

x {4 m 2 G 2 E {q 2 ) - q 2 G 2 M (q 2 )} , (10) 

where G E (t) and Gm(^) are the Sachs form factors. 

The elastic contributions to the time-ordered ampli¬ 
tudes T \, T 2 cannot be specified as easily. In perturba¬ 
tion theory, they are usually referred to as Born terms 
and it is not a trivial matter to specify them at higher 
orders of the calculation. In effective low-energy theo¬ 
ries, the decomposition into a Born term and a ’struc¬ 
ture part’ is not a simple matter, either. For a detailed 
discussion of these aspects, we refer to P%H4U1 . In the 
framework of dispersion theory, however, the decom¬ 
position is unambiguous. The reason is that analytic 
functions are fully determined by their singularities and 
their asymptotic behaviour: dispersion theory provides 
a representation of the amplitudes in terms of its singu¬ 
larities. In our framework, this representation is given 
by the dispersion relations 0 and the sum rule 0. The 
elastic contribution is the part of the amplitude which 
is generated by the singularities due to the elastic in¬ 
termediate states. These are specified in (USD- Accord¬ 
ingly, the elastic parts of Tf, T 2 are obtained by simply 
replacing Vi. V 2 with Vf 1 . Vf 1 and dropping the Regge 
contributions. In the case of X 2 , this leads to 




Inserting the explicit expressions for the elastic contri¬ 
butions to the structure functions, we obtain 


/riel, 2i _ _ J _ 

1 (4m 2 v 2 — q 4 )(4m 2 — q 2 ) 

x {G%:(q 2 ) ~ G 2 M (q 2 )} , 


(4m 2 is 2 — q 4 )(4m 2 — q 2 ) 
x {4m 2 G 2 E (q 2 ) - q 2 G 2 M (q 2 )} . (15) 

Both functions tend to zero when v becomes large: by 
construction, the elastic part of Ti(v, q 2 ) does not con¬ 
tain a singularity at infinity. Moreover, as demonstrated 
in appendix 10 even taken by itself, the elastic contri¬ 
butions can be represented in manifestly causal form. 

The explicit expression for the elastic part of the 
subtraction function, 

Siiq 2 ) = J7J 2 2 \(G 2 E (q 2 ) — G 2 m ( q 2 )) , (16) 

q z (4m z — q z ) 

exclusively involves the form factors, which are known 
very precisely. 

3.2 Inelastic part 

We refer to the remainder as the inelastic part of the 
amplitude: 

T l (v 1 q 2 ) = Tf(v,q 2 )+Tr\v,q 2 ), * = 1,2. (17) 

In contrast to the elastic part, which contains the poles 
generated by the elastic intermediate states and is sin¬ 
gular at the origin, the inelastic part is regular there. 
At high energies, the converse is true: while the elas¬ 
tic part tends to zero, the inelastic part includes the 
contributions from the Reggeons, which are singular at 
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infinity. In particular, the sum rule for the inelastic part 
of the subtraction function reads: 


Expressed in terms of these, S'™ el (g 2 ) and S™ el (q 2 ) are 
given by 


S™\q 2 )=T 1 R (0,q 2 ) 

+2 r± {Vl ( Vt t)- v «( V ,f)} , 

J Vt.-h V 


(18) 


where v t h = M n + ( M 2 — q 2 )/2m denotes the inelastic 
threshold. The dispersive representation for the inelas¬ 
tic part of Ti(v,q 2 ) then becomes 


Tj nel (z/, q 2 ) = S[ ne \q 2 


+ 2v 2 



dv' 

— 


biK,g 2 ) 

yll _ yl _ jg 


(19) 


In the case of T 2 (v, q 2 ), a subtraction is not needed. The 
contribution from the elastic intermediate state to the 
dispersion integral in (O coincides with the expression 
for T^iy, q 2 ) in (flbl) . Removing this part, which is even 
more singular at the origin than T®*(i/, g 2 ), we obtain 
the following representation for the inelastic part: 


ninel 


(v,Q 2 ) 



V n _ yl _ jg 


( 20 ) 


3.3 Subtraction function in terms of cross sections 


Sr\q 2 ) = — £i(Q 2 ), 

®em 

S™\q 2 ) =—£2(Q 2 ), 

^em 

r 1 (Q 2 ) = -r T (Q 2 ) + r 1 i (Q 2 ), 

£ 2 (Q 2 ) = £ T (Q 2 ) + ^(Q 2 ). (26) 

While the transverse parts of £i(Q 2 ) and £ 2 {Q 2 ) only 
differ in sign, the longitudinal parts are quite differ¬ 
ent. Regge asymptotics implies that <tt as well as a l 
grow in proportion to v a ~ x . Accordingly, the integral 
£ T (Q 2 ) converges - it represents a generalization of 
the integral relevant for the Baldin sum rule to Q 2 ^ 0 
(cf. subsection l4.2l) . While £ 2 (Q 2 ) is dominated by the 
contributions from the low-energy region and rapidly 
converges as well, it is essential that Reggeon exchange 
be accounted for in £^(Q 2 ). 

We are assuming that, at high energies, the longi¬ 
tudinal cross section can be approximated with a rep¬ 
resentation of the form 


cr£(iz, Q ) = 8n 2 a e 


2 mv — Q 2 


x ^2 /3 a (Q 2 )(2mv - Q 2 +m 2 

a>0 


S 0 ) a • (27) 


The structure function Vi(^, q 2 ) is a linear combination 
of the transverse and longitudinal cross sections, see 
appendix [A] 


Vi(v,q 2 ) = k(v,Q 2 ){a L (v,Q 2 ) - a T {v,Q 2 )} , 


Gl{v,Q ) = Ty 2 (T L{y i Q ), 


Q 2 


k(v, Q 2 ) 


I v — Q 2 /2m 
2tt 2 v(v 2 + Q 2 ) 


( 21 ) 


The representation of the subtraction function thus in¬ 
volves integrals over the transverse and longitudinal 
cross sections. For <S'i nel (<Z 2 ) and S^ el (q 2 ), the follow¬ 
ing integrals are relevant: 


At Q 2 = 0, a Reggeon term proportional to v a in V\ 
corresponds to a contribution to dx that is proportional 
to v a+1 . For nonzero values of Q 2 , however, the factor 
in front of the sum implies that the corresponding cross 
section contains sub-leading contributions. As discussed 
at the end of section [21 the specific form used for the 
Regge parameterization is not essential - as long as it 
satisfies a once-subtracted dispersion relation and cor¬ 
rectly represents the asymptotic behaviour of the phys¬ 
ical cross section. We stick to the one specified in ©, 
which leads to (EH). 


3.4 Chiral expansion 


z t (Q 2 ) = j 

r-OO 

1 dv k(v, Q 2 ) a T (v : Q 2 ), 

I'th 

(22) 

HQ 2 ) = - 

^Tf( 0,q 2 ) 
m 

nOO 


+ J 

' dv k(v , Q 2 ) Ad L (v, Q 2 ), 

I'th 

/»oo 

(23) 

z* (Q 2 ) = j 

1 dv k(v,Q 2 ) a L (v,Q 2 ), 

^th 

(24) 

l{v, Q 2 ) = a L (v, Q 2 ) - erf (v, Q 2 ). 

(25) 


Chiral perturbation theory (yPT) exploits the fact that 
in the limit m u ,md —> 0 (at fixed Aqcd, m 8 ,..., nit) 
QCD acquires an exact chiral symmetry, which strongly 
constrains the low-energy properties of the amplitudes. 
The chiral perturbation series provides a representation 
of the quantities of interest in powers of momenta and 
quark masses. In the chiral limit, the pion is a massless 
particle, but when the quark masses m u , m<i are turned 
on, the pion picks up mass in proportion to the square 
root thereof, M 2 = (m u + rrid)B + 0{m 2 logm g ). 
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In the context of the present paper, we only need the 
chiral expansion of the form factors Ge(<1 2 ), Gm(q 2 ), 
and of the functions Si(q 2 ), S 2 (q 2 )- These quantities 
involve a single momentum variable, q 2 . As we work 
in the isospin limit m u = rrid , the corresponding chi¬ 
ral perturbation series involves an expansion in the two 
variables M n and q 2 . The series can be ordered in pow¬ 
ers of M n \ the coefficients then depend on the ratio 


AM 2 ’ 


(28) 


which counts as a quantity of 0(1)- In contrast to the 
straightforward Taylor series in powers of q 2 , the chiral 
expansion is able to cope with the infrared singularities 
generated by the pions. 

To leading order in the chiral expansion, the infrared 
singularities are described by a set of one-loop graphs 
of the effective theory [41] . In the case of the magnetic 
Sachs form factor, for instance, the evaluation of the 
relevant graphs within Heavy Baryon yPT leads to the 
following expression for the first non-leading term in 
the chiral expansion [42] FI 


G P M (q 2 )=V P 


g\mMv 

16ttF 2 


(1 + t) 


arctan yj r 

Vt 


gUq 2 ) 


+ 0{M 2 log M n ), 

n : g 2 A mM n fn , ^ arctan Vr 

+ Tfa?ri (1+T) ^^ 


+ 0(M 2 log M n ) ■ 



Up to and including 0(M n ), the magnetic form fac¬ 
tor can thus be represented in terms of the mag¬ 
netic moment /i, the pion decay constant, F v = 
92.21(14) MeV [H], and the nucleon matrix element of 
the axial charge, gA = 1.2723(23) 23- The formula 
shows that, up to higher-order contributions, the sin¬ 
gularity is described by a function of the ratio r = 
(—q 2 )/AM 2 : the scale is set by the pion mass, not by 
Aqcd- The presence of a scale that disappears in the 
chiral limit also manifests itself in the slope of the form 
factor at q 2 = 0, i.e. in the magnetic radius: the above 
representation shows that the chiral expansion of the 
magnetic radii of proton and neutron starts with a term 
of 0(l/M n ). 

5 Note that the range of validity of the representation l[29[) is 
limited. The mr intermediate states generate a branch point 
in the form factors at q 2 = AM 2 , which corresponds to r = 
— 1. While relativistic formulations of Baryon xPT do cover 
this region, an infinite series of Heavy Baryon \PT graphs 
contributes in the vicinity of that point, more precisely in the 
region where r+ 1 is small, of 0(M 2 /m 2 ) 1431 . In the present 
paper, however, we make use of the chiral expansion only near 
r = 0, where the nonrelativistic framework is adequate. 


The low-energy behaviour of the electric Sachs form 
factors is less singular: 

G p E (q 2 ) = 1 + 0(M 2 logM w ), 

G n E (q 2 ) = 0(M 2 log M n ). (30) 


Accordingly, the chiral expansion of the electric radii 
does not start with a term of 0(1/M n ), but with a 
chiral logarithm, comparable to the situation with the 
charge radius of the pion. 

The subtraction function also diverges if the chiral 
limit is taken at a fixed value of the ratio q 2 /M 2 : the 
leading term in the chiral expansion of S'( nel (q 2 ) is of 
order 1/M V and is determined by F n and gA as well [46] : 


sr\q 2 ) 


g 2 A m 

64-7T F 2 M^t 


+ 0(log M n ). 



arctan y/r 1 

Vt j 


(31) 


The expansion of the analogous term in X 2 starts 

with [46] : 


^'inel/ , ^,2 


(g 2 ) = - 


g 2 A m 


64-7T F%MirT 
+ 0(log Mir) ■ 


1 - (1 + At) 


arctan yfr 

Vt 


(32) 


In either case, the leading term is the same for proton 
and neutron - for Sf' el (q 2 ) and S?f el (q 2 ), the chiral ex¬ 
pansion of the difference between proton and neutron 
only starts at OilogM^). 


4 Nucleon polarizabilities 


4.1 Low-energy theorems 


In contrast to the elastic parts, which are singular at the 
origin, the inelastic contributions to Ti(i/,q 2 ),T 2 (v,q 2 ) 
do admit a Taylor series expansion in powers of v and 
q 2 . Two low-energy theorems relate the leading terms 
in this expansion to the polarizabilities of the nucleon. 
The theorems amount to rather nontrivial statements, 
because the functions T\(y, q 2 ), T 2 {iy, q 2 ) represent the 
virtual Compton scattering amplitude in the forward 
direction, while the experimental determination of the 
polarizabilities relies on real Compton scattering at 
nonzero scattering angle. A concise derivation is given 
in appendix [B] 

In the above notation, the low-energy theorems take 
the simple form: 


S[ ne \ 0) = - 


4m 2 


■ Pm , 


sr\ 0 ) =- (ole + Pm), 

^em 


(33) 

(34) 
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OL E 


a E 

0M 


0M 

ole + Pm 


p 10.65(0.50) 

m 

0.55 

3.15(0.50) 

m 

-0.55 

13.8(4) 

m 

n 11.55(1.50) 

m 

0.62 

3.65(1.50) 

m 

-0.62 

15.2(4) 

ssa 

p — n —0.9(1.6) 


-0.08 

—0.5(1.6) 


0.08 

-1.4(6) 



Table 1 Experimental values of the nucleon polarizabilities, in units of 10 4 fm 3 , as determined from EFT extractions in 
Compton scattering I47II49II5TI and analyses of the Baldin sum rule [481150] (see also |371l52p . The latter results were imposed 
in mm, so that the quoted errors for ole and Pm are anticorrelated. 


where n is the anomalous magnetic moment, qb and 
Pm are the electric and magnetic polarizabilities of the 
particle, and a e m is the fine structure constant. These 
relations show that the polarizabilities contain an elas¬ 
tic as well as an inelastic part, while their sum, cxe+Pm, 
is purely inelastic: 


4m 3 


P M = - 


4 TO 3 


(35) 


Table [Q shows that the elastic parts only represent a 
small fraction of the polarizabilities. 


4.2 Sum rules for the polarizabilities 


The low-energy theorem (1531) thus represents the famil¬ 
iar Baldin sum rule [54] . The integral occurring here 
is a limiting case of the quantity S T (Q 2 ) introduced 
in (1221) : the Baldin sum rule amounts to 

a E +p M = ^ T (0). (39) 


Comparison with (1371) shows that the electric polariz¬ 
ability obeys a sum rule that exclusively involves the 
longitudinal cross section and the anomalous magnetic 
moment: 


ole — (0) + 


4m 3 


(40) 


The left hand side of the low-energy theorem (l33l) rep¬ 
resents the inelastic part of the subtraction function at 
q 2 = 0: 


a, 


PZ el = -— Si (0). 

m 


(36) 


The representation for the subtraction function in (USD 
thus amounts to a sum rule for the inelastic part of 
the magnetic polarizability. Adding the elastic contri¬ 
bution, the sum rule takes the form 

Pm = U T ( 0 )-^( 0 )-^- ( 37 ) 

To our knowledge, this sum rule is new. It states that, 
in the absence of fixed poles, the magnetic polarizabili¬ 
ties of proton and neutron are determined by the cross 
sections for photo- and electroproduction. If the ampli¬ 
tude Ti(v,q 2 ) were to obey an unsubtracted dispersion 
relation, the Regge terms in the expression for U^(Q 2 ) 
could be dropped, so that the sum rule would reduce 
to the one proposed in [53]. Regge asymptotics implies 
that a subtraction is needed, but if the Reggeon trajec¬ 
tories and residues are known, the subtraction can be 
expressed in terms of these. 

Evaluating the dispersive representation (EOD at v = 
q 2 = 0, we obtain 


fO O 7 

S'r\ 0 ) = 2 / —v 2 (v,0) 


m 

27r 2 a em 



dv 


ctot iy) ■ 


(38) 


5 Numerical analysis 

5.1 Experimental information 

We evaluate the cross section integrals on the following 
basis H 

W < 1.3 . At low energies, the resonance Z\(1232) 
generates the most important inelastic contribution. It 
decays almost exclusively into ttN final states which 
have been thoroughly explored. The SAID, MAID, 
Dubna-Mainz-Taipei (DMT), and chiral-MAID col¬ 
laborations provide pion photo- and electroproduction 
cross sections into these channels [55IRT31 FI For W < 1.3 
and real photons ( Q 2 = 0), the transverse cross section 
is well approximated by the sum over these contribu¬ 
tions. In particular, the representations we are using 
are consistent with isospin symmetry, which implies 
that the contributions from the A to the proton and 
neutron cross sections are the same up to symmetry¬ 
breaking effects of 0(m u — rad, Q e m) which are expected 
to be very small. Moreover, as seen from figure [1] (left 
panel), the A dominates in the transverse cross sec¬ 
tions and gives very small contributions to the longi¬ 
tudinal ones. This property is directly related to the 

6 Throughout, the numerical values of W and Q refer to GeV 
units, the cross sections are given in fib , while the polariz¬ 
abilities as well as the cross section integrals S are expressed 
in units of 10 -4 fm 3 . 

7 The photoproduction cross sections are also provided by the 
Bonn-Gatchina collaboration | 64ll65[ . 
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Fig. 1 Cross sections ctt and (tl = v 2 /Q 2 a l in the Zl-region, for Q 2 = 0. 


smallness of the C 2 Coulomb quadrupole form factor 
for the AN 7 * transition. In the non-relativistic quark 
model, where both the nucleon and the A are zero- 
orbital-momentum three-quark states, this form factor, 
as well as the one of the E2 electric quadrupole, vanish 
altogether pTHlRTSl . 

The comparison of the full lines in the two panels of 
figure [T] shows that, in the region where the A generates 
the dominant contribution, the transverse cross sections 
for proton and neutron are indeed nearly the same: the 
differences are smaller than the individual terms by an 
entire order of magnitude [ 66 ] . For the polarizabilities, 
the behaviour of the ratio ctl/Q 2 in the limit Q 2 —> 0 is 
relevant. Since MAID and DMT offer a representation 
also for this quantity, these parameterizations are par¬ 
ticularly convenient for us. For definiteness, we identify 
the central values of the cross sections in the region 
W < 1.3 with the average of MAID and DMT, ab¬ 
breviated as MD: omd = ^(umaid + odmt)- As far as 
the proton cross sections are concerned, the results ob¬ 
tained with SAID, MAID, and DMT are practically the 
same, but figure [ 2 ] shows that for the small differences 
between proton and neutron, this is not the case. The 
uncertainties in the input used for the cross sections 
do affect our numerical results and will be discussed 
together with these. 

1.3 < W < 3. In the intermediate region, we rely 
on the work of Bosted and Christy (BC), who pro¬ 
vide parameterizations of the transverse and longitudi¬ 
nal proton and neutron cross sections in the resonance 
region, m + M 7 r < W < 3.2, in the range 0 < Q 2 < 
8 [691170] . These contain a wealth of information, but 
suffer from a number of shortcomings. In particular, 
their fit to the data is carried out under the assumption 
that the ratio ul/pt is the same for proton and neu¬ 
tron. An experimental analysis that does not rely on 


this assumption would be most welcome. Second, the 
parameterization does not properly cover the region of 
very small photon virtualities (cf. [71]): (a) The alge¬ 
braic form of the representation used for ax implies that 
the quantity dy, = (Jlv 2 /Q 2 disappears when Q tends 
to zero instead of approaching a nonzero limiting value, 
(b) Isospin symmetry implies that the contributions of 
the resonance Z\(1232) to proton and neutron are the 
same, but, as noted in [36] . the BC-parameterization 
does not respect this symmetry to the expected ac¬ 
curacy. (c) The parameterization of the contribution 
from the resonance IV(1530) exhibits an unphysical de¬ 
pendence on Q 2 : in the tiny interval 0 < Q 2 < 0.001, 
the contribution from this resonance to the transverse 
cross section of the proton varies by about 40%. Al¬ 
though this artefact only manifests itself at very small 
values of Q 2 , it seriously affects our calculation because 
the results obtained for the polarizabilities depend on 
whether we simply evaluate the sum rules at Q 2 = 0 or 
use very small positive values of Q 2 - for the physical 
cross sections, a difference of this sort cannot arise. 

In the interval 1.3 < W < 3, we use the follow¬ 
ing crude estimate for mean values and errors: (i) The 
central value is identified with the result obtained with 
the BC-parameterization. (ii) In order to wash out the 
spikes occurring at very small values of Q 2 , we assign an 
8% uncertainty to the BC-representation of the proton 
cross sections: Aa v = 0.08 <j p . (iii) Since the difference 
between the proton and neutron cross sections is much 
smaller than the individual terms, small relative errors 
in the latter can generate large relative errors in the 
difference. For this reason, we use the same error es¬ 
timate for <7 p ~ n as for the individual terms, i.e. work 
with Aa p ~ n = 0.08 <j p . 

The comparison of the representations for the dif¬ 
ference between the proton and neutron cross sections 
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Fig. 2 Consistency check at the transition point W = 1.3. The plot compares the representations of MAID and DMT used 
below that point with the BC-parametrization used above it. The difference between MAID and DMT and the band attached 
to BC represent an estimate of the uncertainties to be attached to these parameterizations. As discussed in the text, the picture 
implies that these parameterizations provide a coherent framework only for Q 2 > 0.5. 


used below and above W = 1.3 offers a consistency test 
on our calculations. Figure [2] compares the representa¬ 
tions of MAID and DMT with the uncertainty band 
attached to BC at the transition point. The left panel 
shows that the representations for the difference of the 
transverse cross sections used below and above that 
point agree with one another only for Q 2 > 0.5. The 
problem arises from the deficiencies mentioned above, 
which prevent us from reliably evaluating the cross sec¬ 
tion integrals at low values of Q 2 . The right panel shows 
that the uncertainties in the difference of the longitu¬ 
dinal cross sections are considerable, but within these, 
the representations used are coherent. 

W > 3. We estimate the contributions from higher 
energies with the representation of Alwall and Ingcl- 
man m • It is based on the vector-meson-dominance 
model of Sakurai, Schildknecht, Donnachie, and Land- 
shoff [7M75] and offers a parameterization of the trans¬ 
verse and longitudinal cross sections of the form 

= p £(Q V p_1 + , 

= pp(Q 2 )s ap ~ 1 + Pp(Q 2 )s olR ~ 1 , (41) 

where s = W 2 is the square of the center-of-mass en¬ 
ergy. The Pomeron cut is approximated by a Regge pole 
at ap = 1.091, while the Reggeons with the quantum 
numbers of / and <22 are lumped together in a single 
contribution with ap = 0.55. 

The Pomeron residues of proton and neutron are 
the same: 

Pp(Q 2 r = pp(Q 2 r, Pp(Q 2 r = Pp(q 2 ) p ■ m 

For the remainder, we follow [3], invoke 517(3), and 
stick to the value of the D/F ratio quoted there (for 


the definition of the Regge couplings D and F and a 
review of their determination, we refer to El): 


Pr(Q 2 Y 

Z 


ZP t r {Q 2 Y, Pr(Q 2 ) 


6F-4D 
9 F-D 


0.74. 


ZPk(Q*Y, 

(43) 


The parameterizations for the structure function F 2 
of Capella et al. USED] and for the ratio ctl/ot of 
Sibirtsev et al. EDI provide an alternative Regge repre¬ 
sentation of the cross sections, which we refer to as CS. 
In figure[3j the consistency check made at the transition 
point W = 1.3 is repeated for W = 3. The plot shows 
that, for Q 2 < 1.4, the central representation of Al is 
indeed contained in the uncertainty band attached to 
BC, while the one of CS runs above it. The comparison 
indicates that, at low values of Q 2 , working with Al 
yields a coherent picture, while with CS this is not the 
case. For Q 2 > 2, however, the situation is reversed: 
there, the Al-representation yields values for the differ¬ 
ence between the transverse cross sections that are too 
small while the CS-representation is consistent with the 
values obtained from BC. This confirms the conclusion 
reached in m-- the above form of the Al-representation 
applies as it stands only for Q 2 < 1. At higher values of 
Q 2 1 the parameterization underestimates the size of the 
structure function F 2 and further contributions have to 
be added for the vector-meson-dominance formulae to 
become compatible with the observed behaviour. Since 
we do not account for these and the uncertainties we 
attach to the central representation do not cover the 
gap, the input we are working with becomes incoher¬ 
ent for Q 2 > 2. The right panel, on the other hand, 
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Fig. 3 Consistency check at the transition point W = 3, see main text. 


shows that the representations we are using for the lon¬ 
gitudinal cross section do survive the consistency test, 
irrespective of the value of Q 2 . 

Note that we are discussing the properties of the dif¬ 
ference between the proton and neutron cross sections. 
The main problem here is that all of the well-established 
properties of the proton cross sections drop out when 
taking the difference between proton and neutron. High 
precision is required to measure the remainder, in par¬ 
ticular also at high energies, where the Pomeron dom¬ 
inates the scenery. Also, since the longitudinal cross 
section is significantly smaller than the transverse one, 
pinning it down accurately is notoriously difficult. In 
both of the above representations, the ratio is 

taken to be energy-independent^ This appears to be 
consistent with experiment, but since we are not aware 
of a theoretical explanation, a test of the validity of 
this assumption would be very useful. For recent appli¬ 
cations of these representations to the amplitudes under 
consideration we refer to [7TU8T1I82] . 


5.2 Evaluation of £ T and 27 2 for the proton 

We start the discussion of the cross section integrals 
with the one over the transverse cross section, £ T (Q 2 ), 
which is specified in (E^l) . The value at the origin is 
relevant for the sum of the electric and magnetic po¬ 
larizabilities, A’ T (0) = otE + /3m- Since the longitu¬ 
dinal cross section vanishes at Q 2 = 0, the function 
A 2 (Q 2 ) = £ T (Q 2 ) + 2/f(Q 2 ) takes the same value 
there. In fact, figure [4] shows that 27 2 (Q 2 ) is dominated 

8 For the representation ED- this property implies dp/dp = 
dp //3p- In the notation of [72], it corresponds to By / Ay = 
B-y/A-y, for V = 


by the transverse part also at nonzero virtuality - the 
longitudinal part amounts to a modest correction. 

As pointed out in EL], the structure function 
Fi(x, Q 2 ) can also be used to continue the integral rel¬ 
evant for the Baldin sum rule to nonzero values of Q 2 : 

ZfAQ 2 ) = r'dxxF, 

V Jo 

1 r °° rl 

= ^l th i("-Q 2 / 2"W- (44) 

Since the integrand differs from the one relevant for 
£ T (Q 2 ) only by the factor l + Q 2 /is 2 , the quantity £p x 
also reduces to qe + /3m when Q 2 vanishes, but drops 
off somewhat less rapidly when Q 2 grows. 

The lines for £ T , 27 2 , and £p 1 in figure [4] are ob¬ 
tained by using the parameterizations specified in sub¬ 
section I5TT1 As stated there, the contributions from the 
region W < 1.3 are evaluated with the mean of MAID 
and DMT, but we could just as well have used SAID 
on this plot, the difference would barely be visible. 

In [?T], the function £pAQ 2 ) is instead evaluated 
with the BC-parameterization, also in the region of the 
Z\-resonance. This leads to the behaviour indicated by 
the dash-dotted line labelled £p 1 . The topmost line, 
which is obtained by evaluating the same formula with 
the MD-parameterization, is higher by about 0.8 units. 
The difference is closely related to the fact that the BC- 
parameterization does not respect isospin symmetry to 
the expected accuracy (see the discussion in subsec¬ 
tion 15.81) . 

As pointed out by Bernard, Kaiser, and 
Meifiner [83], xPT neatly explains the size of the 
combination of polarizabilities occurring in the Baldin 
sum rule. The parameter-free expression (1321) for the 
leading term in the chiral perturbation series of A 2 ((5 2 ) 
is shown as a dashed line. The comparison with the 
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Fig. 4 Cross section integrals related to the Baldin sum rule 
for the proton (numerical values in units of 10 -4 fm 3 ). At 
Q 2 = 0, all of the quantities shown reduce to oP E + f) p M . The 
plot focuses on small values of Q 2 , where the behaviour is 
dominated by the Nambu-Goldstone bosons - in the chiral 
limit these generate an infrared singularity. The short-dashed 
line shows the parameter-free result obtained from xPT a t 
leading order. The cross section integrals are evaluated with 
the parameterization specified in subsection 15.11 except for 
Sp 1 , where the contribution from the region of the Delta is 
calculated with BC instead of MD. 


experimental result for + 0m shows that, at small 
values of Q 2 , the leading term of the chiral series 
dominates. In the limit Q 2 —> 0, this term reduces to 


ole + 0m 


1 

192ttF%Mt T ' 


(45) 


In the chiral limit this formula diverges in inverse 
proportion to ill,: if the quarks are taken massless, 
Tr>> 9 2 ) contains an infrared singularity at v = q 2 = 

0 . 

The same singularity also shows up in the in¬ 
dependence, which exhibits the presence of an unusu¬ 
ally small scale: at leading order of the chiral expansion, 
the function ^(Q 2 ) depends on Q 2 only via the vari¬ 
able r = Q 2 /4M 2 . Hence the scale is set by 2M n rather 
than M p . Figure 0 ] shows that, in reality, Ei(Q 2 ) drops 
even more rapidly, partly on account of the second-sheet 
pole associated with the A, partly due to other higher- 
order contributions of the chiral series snraiHSj. 

The spike seen in figure 0] at tiny values of Q 2 illus¬ 
trates the artefact mentioned in subsection 15.11 which 
concerns the contribution from the resonance 7V(1530): 
if the numerical values of the integrals in the region 
0.002 < Q 2 < 0.005 are fit with a low order polynomial, 
the extrapolation to Q 2 = 0 is higher than the result 
of the direct evaluation at Q 2 = 0, by about 0.4 units. 
Since the experimental information from real Compton 
scattering and from photoproduction is more stringent 



Fig. 5 Cross section integrals relevant for the difference be¬ 
tween proton and neutron. 


than the one from electron scattering, which for these 
very small values of Q 2 necessarily involves extrapola¬ 
tions, we think that the results obtained by evaluating 
the integral over the transverse cross section at Q 2 = 0 
are more reliable. The value obtained there with MAID 
or DMT is (a_E + 0m) p = 14.1, while SAID yields a 
result that is lower by about 0.1 units. The numbers 
obtained at Q 2 = 0 with the parameterizations we are 
using thus agree with the result (a_E + 0m) p = 13.8(4) 
quoted in the review EH, which stems from [48]. 


5.3 E T and E 2 : proton-neutron difference 

Figure [5] shows the difference between the integrals over 
the proton and neutron cross sections. The picture looks 
very different from figured while there, the curves start 
at E cs 14 and rapidly drop with Q 2 , those in figure 0] 
start at E ~ 0 and stay there. Since the integrals under 
consideration are rapidly convergent, the behaviour of 
the cross sections in the resonance region is relevant. 
The reason why not much is left in the difference be¬ 
tween proton and neutron is that, in that region, the 
proton and neutron cross sections are nearly the same. 
In particular, as mentioned in subsection 15.11 isospin 
symmetry implies that the most prominent low-energy 
phenomenon, the A, drops out when taking the differ¬ 
ence between the proton and neutron cross sections. 
The cancellation of the main contributions also mani¬ 
fests itself in the chiral perturbation series: the leading 
terms in E p and E% are large, of order 1 /M n , but the 
coefficients are the same, so that the chiral expansion 
of E p ~ n only starts at 0(1). 

For our cross section integrals to exhibit these fea¬ 
tures, it is essential that the representations we are us- 
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ing in the region of the A respect isospin symmetry. 
The dash-dotted line illustrates the fact that the BC- 
parameterization of the cross sections violates this con¬ 
straint quite strongly: the bump seen around Q 2 ~ 0.1 
arises from the difference between proton and neutron 
which occurs in that parameterization in the region of 
the A. As mentioned above, the difference between the 
parameterizations MD and BC in the region W < 1.3 
also shows up in figure Q] It so happens that the dif¬ 
ference between the results obtained via extrapolation 
from Q 2 > 0.002 and via evaluation at Q 2 = 0 nearly 
cancels the one between the contributions from the re¬ 
gion of the A obtained with BC and with MD, so that 
the number obtained for a p E + fi p M in [81] agrees with 
experiment. 

The spike seen at very small values of Q 2 is about 
twice as large as the one in figure [I] and manifests it¬ 
self much more prominently because the difference be¬ 
tween proton and neutron is an order of magnitude 
smaller than the individual terms. The value obtained 
at Q 2 = 0 is consistent with the experimental result, 

(aE + /3m) p-ti = —1.4(6). 

5.4 Pomeron exchange 

The integrals considered in the preceding two subsec¬ 
tions converge rapidly. Their properties are governed 
by the low-energy behaviour of the cross sections - the 
asymptotic behaviour does not play a significant role. 
For the integral E^(Q 2 ) specified in (1231) . the situa¬ 
tion is very different: for this integral to converge, it is 
essential that the asymptotic behaviour of the longitu¬ 
dinal cross section be known, so that it can properly 
be accounted for. At high energies, the leading contri¬ 
bution stems from Pomeron exchange, which generates 
a branch point at J = 1 in the angular momentum 
plane. In phenomenological parameterizations, such as 
the one specified in m , the branch cut is often ap¬ 
proximated by a Regge pole in the range 1 < ap < 2. 
For this parameterization to have the required asymp¬ 
totic accuracy, it must describe the contribution from 
the Pomeron up to terms that disappear in the limit 
v —> oo. 

The Regge representation we are using to describe 
the asymptotic behaviour of the structure functions 
leads to the parameterization (1271) . In this framework, 
the Pomeron term in (Till) not only generates a leading 
contribution to the cross section with a = ap, but also 
a daughter with a = ap — 1. Furthermore, in contrast 
to the situation with the parameterization of the contri¬ 
butions from the nonleading Reggeons, the value of the 
parameter So does matter here: a change in the value of 
sq generates an asymptotic contribution proportional 


to v ap 1 . If the integral in (E3l) does converge for one 
particular value of So, it diverges for any other value. 

As an illustration of the mathematical problem we 
are facing here, consider a contribution of the form 

AT\{y, q 2 ) = ^(<? 2 )|(si - m 2 - 2 mv - q 2 ) 8 

+ (si — to 2 + 2 mv — q 2 ) <5 |, (46) 

which is free of fixed poles. For v > 0, q 2 < 0, the 
corresponding absorptive part is given by: 

AV^(v,q 2 ) = £(g 2 ) 6>(m 2 + 2 mv + q 2 - Si) 

Z7T 

x (m 2 + 2 mv + q 2 — si) 6 . (47) 

In the limit S —> 0, the modification of the structure 
function disappears, while the change in the time-ordered 
amplitude does not, but takes the form of a fixed-pole 
contribution, ATi(v,q 2 ) —> £(g 2 ), which can have any 
desired value. 

In short: although the hypothesis that the Reggeons 
properly account for the behaviour at large values of 
v uniquely determines the subtraction function even if 
Pomeron exchange contributes, the evaluation of C3l) 
requires knowledge of the asymptotic behaviour to 
an accuracy that is beyond reach. In the absence of 
theoretical information about the properties of the 
Pomeron, we are dealing with what Hadamard [86] 
called an ill-posed problem: in principle, the data do 
determine the solution, but tiny changes in the data 
(structure function) can lead to substantial changes in 
the solution (subtraction function). For this reason, we 
do not discuss the sum rules for the individual polariz¬ 
abilities of proton and neutron any further. 

A model-independent determination of the subtrac¬ 
tion function occurring in the dispersive representa¬ 
tion of the proton Compton amplitude is also of in¬ 
terest in connection with the proton radius puzzle (for 
a recent review see EH)- As pointed out in [88], at 
least part of the discrepancy could be explained if for 
some reason the contribution to the Lamb shift that is 
governed by the virtual Compton scattering amplitude 
were significantly larger than expected. The X'PT analy¬ 
ses pH)ll46ll84l[55] as well as the recent works on effective 
field theory m and finite-energy sum rules E5 were 
largely motivated by this puzzle; an improved knowl¬ 
edge of the subtraction function would be of interest 
also in that context. Unfortunately, however, a major 
breakthrough in the theoretical understanding of the 
Pomeron is required before the sum rule set up above 
could reliably be evaluated. 
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Fig. 6 Cross section integrals relevant for the subtraction 
function. 


5.5 Evaluation of E^ for the proton-neutron difference 

In the present subsection, we focus on the difference 
between proton and neutron, where the Pomeron drops 
out: the asymptotic behaviour of ox(i/, Q 2 ) p ~ n is dom¬ 
inated by the nonleading terms in dm, which grow less 
rapidly with is, so that the problems discussed in the 
preceding subsection do not arise. The analysis of the 
difference is of interest for two reasons: (1) the sum rule 
is obtained under the same premises (absence of fixed 
poles, Reggeon dominance hypothesis) as the Cotting- 
ham formula. Consequently, confronting the result with 
existing experimental information on the polarizabili¬ 
ties, one may test the validity of this hypothesis. (2) 
our result for a p E n is somewhat more accurate than 
the determination based on the current experimental 
information. Combined with the experimental values of 
the polarizabilities of the proton and the Baldin sum 
rule, this yields an improved prediction for the polariz¬ 
abilities of the neutron. 

Figure [G] compares the integrals over the transverse 
and longitudinal cross sections, for the difference be¬ 
tween proton and neutron. The function E T (Q 2 ) p ~ n 
already occurred in figure [3] we are now merely fo¬ 
cusing on a smaller range in the variable Q 2 . The plot 
shows that the integral E^(Q 2 ) p ~ n behaves in a qual¬ 
itatively different way. Both integrals are small, but 
while E t (Q 2 ) exhibits the pronounced spike at Q 2 = 0 
discussed earlier, the dependence on Q 2 of Ei(Q 2 ) is 
dominated by the contribution from the region of the 
A, which is well understood in particular, the MAID 
and DMT representations show nearly the same In¬ 
dependence. Using the mean of the two as central value 
and half of the difference as an estimate for the un¬ 
certainty for the contributions from W < 1.3 would in 


our opinion represent a fair recipe, but to stay on the 
conservative side, we double the error estimate. For the 
value of the integral at Q 2 = 0 this prescription yields 
£('(0)md = —1.4(4). The contributions from interme¬ 
diate energies, 1.3 < W < 3, are small: the estimate 
^’i'(0)bc = 0.2(2) covers the deficiencies of the repre¬ 
sentation used there. Above that range, we use the AI- 
representation, attach an uncertainty of 30% to it, and 
get £f(0) A i = —0.3(1). Adding errors in quadrature, 
we finally obtain 

£f(0) p_n = -1.6(4). (48) 


5.6 Prediction for the polarizabilities of the neutron 

In view of the relation pop . the result (1151) amounts 
to a prediction for the difference between the electric 
polarizabilities of proton and neutron: 

«T n = -1-7(4). (49) 

This is consistent with the current experimental value, 
oP E n = —0.9(1.6), but significantly more precise. The 
numerical result obtained from the Baldin sum rule for 
the difference in the value of ole + Pm between proton 
and neutron, (a# + Pm) p ~ u = —1-4(6), then implies 

P P M n = 0-3(7). (50) 

According to (1551) . this result also determines the value 
of the subtraction function relevant for the self-energy 
difference at Q 2 = 0: 

^inei( 0 )p-n = _ 0 .3(1.2) GeV -2 . (51) 

Finally, combining the current experimental result for 

the electric and magnetic polarizabilities of the proton, 
a p E = 10.65(50) and P^ = 3.15(50), with the numbers 
for a p E ~ n and {aE+PM ) p ~ n , we arrive at a prediction for 
the electric and magnetic polarizabilities of the neutron: 

a£ = 12.3(7), Pm = 2.9(0.9). (52) 

These are also consistent with the current experimental 
values, Om = 11.55(1.5), P E = 3.65(1.50), and more 
precise. 

Note that the procedure used avoids relying on the 
available parameterizations of the transverse cross sec¬ 
tion, which contain sharp spikes at very small values of 
Q 2 that make the evaluation of E T { 0) problematic. We 
make use of the fact that those present in the longitudi¬ 
nal cross section are much milder and allow us to assign 
a meaningful uncertainty to E[( 0). We also emphasize 
that the fluctuations exclusively affect the behaviour 
at small values of Q 2 . For the evaluation of the electro¬ 
magnetic self-energy to be discussed in section [6l these 
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deficiencies are of no concern, because phase space sup¬ 
presses the contributions from the vicinity of the point 

Q 2 = 0. 


5.7 Result for the subtraction function 

According to ( 1551 ) . the inelastic part of the subtrac¬ 
tion function relevant for the self-energy is determined 
by the difference between the integrals £i(Q 2 ) p ~ n and 
£ T {Q 2 ) p ~ n ■ The central values of these integrals are 
shown in figure [6] The narrow band in figure [7] indi¬ 
cates the corresponding result for the subtraction func¬ 
tion. The width of the band is obtained by evaluat¬ 
ing the uncertainties in the contributions arising from 
the three subintervals, separately for the transverse and 
longitudinal contributions, and adding the results in 
quadrature. For better visibility, the vertical axis is 
stretched with the inverse of the dipole form factor, 
N = (1 4- Q 2 /M J) 2 , M 2 = 0.71 GeV 2 . As discussed in 
subsection EU the region Q 2 < 0.5 contains unphysi¬ 
cal fluctuations - this is why we chop the uncertainty 
band off there. Note also that, although the calcula¬ 
tion returns reasonable results even at Q 2 = 2, it is 
not reliable there, because it does not account for the 
contributions by which the Al-parameterization needs 
to be supplemented in order to agree with experiment 
at those values of Q 2 (see subsection 15. II) . 

The figure also indicates the value S’J nel (0) p-ra = 
1.0(2.7) obtained from the current experimental result 
for , as well as our prediction in (IFTl) . These num¬ 
bers concern the value of the subtraction function at 
Q 2 = 0, but are slightly displaced for better visibility. 


5.8 Comparison with previous work 

Recently, Walker-Loud, Carlson, and Miller m pro¬ 
posed a simple ansatz for the subtraction function. In 
our notation, their proposal amounts to 


<S'wcm(9 2 ) = — 


1 


rrin 


2 9 

rag — q z 


m/3 M 


+ -{G 2 M {q 2 )-F 2 D {q 2 )}. 


(53) 


The singularity at q 2 = 0 arises from the elastic con¬ 
tribution in (1151) . The corresponding expression for the 
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Fig. 7 Momentum dependence of the subtraction function 
(GeV units). For better visibility, the vertical axis is stretched 
with the inverse of the dipole form factor, N = (1 + Q 2 /M 2 ) 2 , 
M 2 = 0.71 GeV 2 . In this normalization, the ansatz proposed 
in [31; (WCM) represents a broad band of nearly constant 
width, determined by the experimental value of the difference 
between the magnetic polarizabilities of proton and neutron. 
The curves are drawn for the current experimental value, 
which is indicated by the error bar on the left and concerns 
the value at Q 2 = 0, but is displaced to make it visible. The 
range obtained with the model in 1361 (ESTY) starts with 
the same width at Q 2 = 0, but shrinks as Q 2 grows. The 
comparatively narrow third band represents our work. We do 
not present an error estimate in the region 0 < Q 2 < 0.5, 
because there, our results are sensitive to the inadequacies of 
the parameterizations used for the cross sections, but we do 
show our prediction for the value of the subtraction function 
at Q 2 = 0. 



inelastic part of the subtraction functional 


cinel 

°WCM 


(q 2 ) 


_ f ml \ 2 mpM 

V TO 0 — q 2 ) «em 

4 m 2 {G E {q 2 ) - G M {q 2 )} 2 
(4m 2 — q 2 ) 2 


(54) 


is regular at q 2 = 0 and one readily checks that the 
ansatz is consistent with the low-energy theorem (1361) . 
It amounts to an extrapolation of that formula to 
nonzero values of q 2 7 controlled by the parameter mo. 
In figure [3 this expression is indicated as a broad band 
of nearly constant width. 

The ansatz ( 1551 ) for the subtraction function gen¬ 
erates a logarithmic divergence in the integral ( 1551 ) for 

9 In [31] . a different terminology is used: there, the first term 
in (1541) is referred to as the inelastic part of the subtraction 
function. This is inadequate, as it amounts to counting the 
polarizability as a purely inelastic quantity. The names given 
to the various terms are not of importance physically, but 
when comparing formulae and numerical values, differences 
in nomenclature must be accounted for. 
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the corresponding contribution to the self-energy dif¬ 
ference. As discussed in subsection EEU the self-energy 
difference indeed diverges logarithmically. The diver¬ 
gence is absorbed in the electromagnetic renormaliza¬ 
tion of m u and which is of order e 2 m u , e 2 rrid■ As 
pointed out by Erben et al. J35], the logarithmic diver¬ 
gence generated by the ansatz (1551) is not proportional 
to the masses of the two lightest quarks and can thus 
not be absorbed in their renormalization: the particu¬ 
lar extrapolation proposed in m is not consistent with 
the short-distance properties of QCD. The variant pro¬ 
posed in [36], 


oinel 

d E sty 


(9 2 ) 


_ / m\-cq 2 \ ( to 2 \ m/3 M 
\m\- q 2 ) \ml - q 2 ) a em 

4 m 2 {G E (q 2 ) - G M (q 2 )} 2 
(4 m 2 — q 2 ) 2 


(55) 


repairs this shortcoming, as it disconnects the be¬ 
haviour at small values of q 2 from the asymptotic be¬ 
haviour. This expression is represented by the central 
band that gradually shrinks if Q 2 increases. 

The explicit choice made in m for the coefficient 
c implicitly assumes that the contribution from the 
subtracted dispersion integral, m^ lsp , stays finite when 
the cut-off is removed, so that the logarithmic diver¬ 
gence then exclusively arises from the term m®, which 
stems from the subtraction function. As discussed in 
detail in [3], however, the deep inelastic region also con¬ 
tributes to the coefficient of the logarithmic divergence. 
The scaling violations do not extinguish this contribu¬ 
tion ESI- Hence the choice made for c cannot be taken 
literally, but it does have the proper quark mass fac¬ 
tors, so that the divergence arising from the subtraction 
function is suppressed. Since the authors cut the inte¬ 
gral over the subtraction function off at A 2 = 2 GeV 2 , 
it barely makes any difference whether c is set equal to 
zero or taken from [33:. In fact, one of the variants of 
the model studied in [35] does correspond to c = 0 . 


6 Self-energy 

6.1 Cottingham formula 

The electromagnetic self-energy of a hadron diverges 
logarithmically. To first order in a em the renormalized 
electromagnetic Lagrangian requires counter terms pro¬ 
portional to the operators 1 , qq , and Oq = G* v G a,iV : 

£em = -y J d 4 yD A (x - y)Tj tl (x)j v (y) + AEt 

+ ^ 99 _ , (56) 

q=u,d,... ^ 


where Da{x) is the regularized photon propagator in 
coordinate space. The counter term proportional to the 
unit operator does not contribute to the self-energy. 
The remainder is determined by the renormalization of 
the quark masses and of the coupling constant g re¬ 
quired by the electromagnetic interaction. To leading 
order, these are given by (see for instance ED): 


3e 2 A 2 

5lJlq= 16 ^ log ^ Q l m i 


Sg = - 


e 2 g 3 


log 


A 2 


2567 r 4 m y 


E 


Q\ 


q—iL,i 


(57) 


The form of the regularization used for the photon prop¬ 
agator is irrelevant - it exclusively affects the value of 
the running scale y. 

The proton and neutron matrix elements of the op¬ 
erator (1551) lead to a version of the Cottingham for¬ 
mula m that is valid in QCD: 

= 2 m( 27 r) 4 / d4 ^A(g 2 ){3q 2 T 1 + {2v 2 + q 2 )T 2 } 

+ counter terms . (58) 


It represents the electromagnetic self-energy in terms 
of the time-ordered amplitudes Tf and T 2 specified in 
appendix [A] 


6.2 Elastic part of the self-energy 

Analogously to the electric and magnetic polarizabili¬ 
ties, the self-energy also consists of an elastic and an 
inelastic part, 

m 7 = m ® 1 + m“ el . (59) 

The contribution from the elastic intermediate 
states remains finite when the cut-off is removed. It 
is obtained by replacing T\ : T 2 with the elastic parts 
Tf 1 , T r f. which are given explicitly in (1T51) . and replac¬ 
ing Dyi{q 2 ) with the full photon propagator, D(q 2 ) = 
{—q 2 — *e) _1 . With a Wick rotation, the expression can 
be brought to the form 

nO O 

< = ^3 j dQ 2 Q 2 {fi vf(-Q 2 ) + f 2 vf(-Q 2 )} , 

f 2 = (l-2y)jl+±-+2y, (60) 

where vf 1 ^ 2 ) and vf(q 2 ) represent the sums of squares 
of form factors specified in (flOl) . The variable y stands 
for y = v 2 IQ 2 . For the elastic contribution, which is 
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concentrated to the line Q 2 = 2mv , we have y = 
Q 2 /Am 2 . In [J], the dipole approximation for the Sachs 
form factors was used, which yields to® 1 = 0.63 MeV 
for the proton and —0.13 MeV for the neutron, so that 
the elastic contribution to the self-energy difference 
amounts to (m^) p ~ n = 0.76 MeV. 

In the meantime, the precision to which the form 
factors are known has increased significantly. For a thor¬ 
ough review of the experimental information, we refer 
to m- The above estimates of the elastic contributions 
to the proton and neutron self-energies do receive sig¬ 
nificant corrections, but the difference between proton 
and neutron is affected by less than 0.02 MeV. Com¬ 
pared to the uncertainties in the contributions arising 
from the deep inelastic region, the departures from the 
dipole approximation are too small to matter. 


6.3 Inelastic part of the self-energy 

The inelastic part receives three distinct contributions: 

m”® 1 = TO^ + mf p + TO®*. (61) 

The term arises from the subtraction function 
gm® 1 ^ 2 ), m disp gj ven by a dispersion integral over 
the structure functions, and to®* accounts for the fact 
that the electromagnetic interaction renormalizes the 
quark masses as well as the coupling constant of QCD. 
In the above discussion of the polarizabilities, renormal¬ 
ization did not play any role, because these concern the 
properties of Ti, T 2 at low energies. In fact, the inelastic 
part of the magnetic polarizability exclusively picks up 
the contribution from the subtraction function speci¬ 
fied in (1M1) . In the decomposition used in (1611) , we have 

/3 inel = j gS ) yjdisp _ pet _ Q 

The term to ^ is obtained by replacing Ti(i/,q 2 ) 
in (l58l) by the subtraction function S'“ el (g 2 ), perform¬ 
ing a Wick rotation, and averaging over the directions 
of the Euclidean momentum. The result reads 


tion for T( nel . The explicit expression reads 


m disp _ “em 


2<n2 


2^xm 


dQ z Q 


dv v 


h~f y \ Vi(v,-Q 2 ) + h V 2 (y, -Q 2 ) } ■ 


(63) 


The term with 3/2 y makes the difference between the 
unsubtracted and subtracted dispersion integral over 
V\\ it removes the leading term in the behaviour of f\ 
when Q 2 is held fixed and v tends to 00 , so that the 
integral over v converges, despite the growth of Vi gen¬ 
erated by Reggeon exchange. On the other hand, when 
Q 2 becomes large, the behaviour in the deep inelastic 
region is relevant. In QCD, the contributions from that 
region diverge logarithmically if the cut-off is removed. 
In (RTT?1) , we have simply cut the integral off at Q 2 = A 2 
- this amounts to a regularization of the photon propa¬ 
gator in Euclidean space: Da(—Q 2 ) = 0{A 2 — Q 2 )/Q 2 . 

In the normalization of the states (IA.2I) . the mass 
shift generated by the counter terms in (1561) is given by 




Jg_ 

Amg 3 


{p\O a \p) • 


(64) 


Neglecting second-order isospin-breaking effects pro¬ 
portional to e 2 (m u — md ), the proton and neutron ma¬ 
trix elements of operators with isospin zero are the 
same. Hence the operators O g , ss, cc, ... drop out in 
the self-energy difference. Moreover, isospin symmetry 
relates the neutron matrix elements of the light quarks 
to those for the proton, e.g. ( k\uu\k) n = (k\dd\k) p . Us¬ 
ing these properties, the contribution from the electro¬ 
magnetic renormalization of the quark masses to the 
self-energy difference can be brought to the form 


CV /\ — 

(m®*) p-rl = - cm (4 to m - m d ) log —(p\uu - dd\p). 

' ZAnm 

(65) 


m * = l^ Jo ^ ^ ^ nel(_g2) ■ (62) 

This term measures the size of the self-energy arising 
from the subtraction function (more precisely, the in¬ 
elastic part thereof - the remainder is included in dm® 1 ). 

The second term on the right of m is obtained 
by replacing the amplitudes Tf, T 2 with their inelastic 
parts T( nel , T™® 1 and dropping the contribution from 
the subtraction function in the dispersive representa¬ 


The formula shows that the coefficient of the logarith¬ 
mic divergence is proportional to the masses of the two 
lightest quarks. In the chiral limit the divergence dis¬ 
appears altogether: if u and d are taken massless, the 
self-energy difference approaches a finite limit if the cut¬ 
off is removed. In reality, the contributions from the 
deep inelastic region do generate a logarithmic diver¬ 
gence, albeit with a small coefficient. An update of the 
analysis performed in [1] is needed to account for the 
scaling violations in the corresponding contributions to 
the renormalized self-energy difference. 
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6.4 Numerical evaluation 

In [5T] , the contribution from the subtraction function 
to the self-energy difference is evaluated with A 2 = 
2GeV 2 . According to ( 1551 ) . the inelastic part of the 
subtraction function is given by the difference between 
two cross section integrals. The part which involves the 
transverse cross sections, E T (Q 2 ), generates a conver¬ 
gent contribution to the formula (l62l) for the self-energy. 
As discussed above, our numerical representation of 
E T (Q 2 ) becomes incoherent at values of Q 2 below 0.5, 
but phase space suppresses that region, so that our es¬ 
timate, m 7 (A T ) ~ — 0.14MeV, should be close to the 
truth (actually, with a coherent representation of the 
available experimental information, this part could be 
evaluated rather accurately, even without cutting the 
integral off). The corresponding integral over the lon¬ 
gitudinal cross section, E^(Q 2 ), is less sensitive to the 
shortcomings of the representation we are using (this is 
why we were able to obtain a rather accurate predic¬ 
tion for the difference between the electric polarizabili¬ 
ties of proton and neutron). Numerically, the contribu¬ 
tion from that integral to the self-energy difference is 
tiny: m^(Z'f') ~ — 0.03MeV. In other words: the con¬ 
tributions from the Reggeons do require a subtraction, 
but taken together with those arising from low energies, 
the entire contribution from the longitudinal cross sec¬ 
tion to the subtraction function generates a negligibly 
small part of the self-energy difference. Together with 
the number for the contributions from the transverse 
cross section given above, we obtain 

TO s = -0.17 MeV. (66) 

This is to be compared with the number obtained by 
instead inserting the expression (1541) in formula (1621) . 
With the central value /3^ n = — 1 used as input in pH] , 
we obtain (m 7 ) WCM = 0.50 MeV. Keeping all other 
parts of the calculation in m as they are, but replacing 
the ansatz for the subtraction function made there with 
our prediction, the numerical result for the self-energy 
difference, m^ CM = 1.30 MeV, is lowered by 0.67 MeV, 
so that the central value becomes m 7 = 0.63 MeV. 
Repeating the exercise with the model of [5B], i.e. re¬ 
placing the expression (1541) by (l55l) . we instead obtain 
(m 7 ) ESTY = 0.20 MeV, so that in this case, the cen¬ 
tral value m ESTY = 1.04 MeV is lowered by 0.37 MeV, 
which leads to m 7 = 0.67 MeV. In either case, the early 
estimate obtained in [J, m 7 L = 0.76(30) is confirmed. 
Comparing their parameterization with recent lattice 
data on the electromagnetic self-energy difference, the 
authors of m and 23 obtain results for the difference 
of the magnetic polarizabilities, 0 p 1 ^ n = —0.87(85) and 
Pai™ = —1-12(40), respectively, which is lower than our 


prediction in (1501) . The difference reflects the fact that, 
in figure [TJ the bands that correspond to their mod¬ 
els run above ours. While these extractions involve a 
model dependence which is difficult to quantify, there 
has recently been progress in the direct calculation of 
the polarizability from the lattice, see [53]. 

Note that the momentum dependence of the sub¬ 
traction function must match the behaviour in the deep 
inelastic region. Taken by itself, the contribution from 
the subtraction function is very sensitive to the choice of 
the cut-off A. As shown in [4|, the term m^ lsp is equally 
sensitive, but the sum of the two contributions is nearly 
independent of A, because the Cottingham formula only 
contains the very weak logarithmic divergence that is 
related to the electromagnetic renormalization of the 
quark masses m u and m<i- As indicated in (1651) . the 
coefficient of the divergence is proportional to these 
masses and hence very small. Also, it does not come ex¬ 
clusively from the subtraction function. The contribu¬ 
tions to m^ lsp arising from the deep inelastic region con¬ 
tribute to the coefficient of the logarithmic divergence 
as well. These were estimated in [4] on the basis of the 
data available at the time, which did not show any vi¬ 
olations of Bjorken scaling. In the meantime, there has 
been considerable progress in understanding the prop¬ 
erties of the structure functions in the deep inelastic 
region and there is very clear evidence for scaling vio¬ 
lations. For a thorough review of these developments, 
we refer to [54]. A corresponding update of the results 
obtained on the basis of the Cottingham formula would 
be of high interest, also in view of the progress made in 
calculating electromagnetic self-energies on the lattice, 
but this goes beyond the scope of the present paper. 


7 Summary and conclusion 

1. Causality relates the imaginary part of the ampli¬ 
tude for Compton scattering on the nucleon in the 
forward direction to the cross section of the process 
e + N —> e + anything. The relation holds for real 
photons as well as virtual photons of spacelike mo¬ 
mentum, q 2 < 0. The spin-averaged forward scat¬ 
tering amplitude involves two invariants, which we 
denote by Ti(u,q 2 ) and T 2 (v,q 2 ). Their imaginary 
parts are determined by the transverse and longi¬ 
tudinal cross sections of electron scattering, <jt and 

&L- 

2. Regge asymptotics implies that only T 2 (^, q 2 ) obeys 
an unsubtracted fixed-g 2 dispersion relation, while 
the one for T\(y,q 2 ) requires a subtraction, which 
represents the value of the amplitude at v = 0: 
5*1 (<? 2 ) = Ti(0,g 2 ). The dispersive representation of 
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the spin-averaged forward Compton scattering am¬ 
plitude thus consists of two parts: an integral over 
the cross sections <jt , ox and an integral over the 
subtraction function S\. The same also holds for 
the Cottingham formula, which represents the elec¬ 
tromagnetic self-energy of the nucleon in terms of 
the spin-averaged forward Compton amplitude. 

3. It had been pointed out long ago [3] that - unless the 
Compton amplitude contains a fixed pole at J = 0 
- the subtraction function is unambiguously deter¬ 
mined by the cross sections of electron scattering. 
We do not know of a proof that the Compton am¬ 
plitude of QCD is free of fixed poles, but assume 
that this is the case and refer to this assumption as 
Reggeon dominance. As briefly discussed in subsec¬ 
tion [OJ the validity of this hypothesis is questioned 
in the literature. Indeed, an analysis of the Comp¬ 
ton amplitude based on first principles that would 
determine the behaviour in the Regge region (high 
energies, low photon virtualities) is not available. If 
the hypothesis were to fail, this would be most inter¬ 
esting, as it would imply that the known contribu¬ 
tions generated by the short-distance singularities 
and the exchange of Reggeons do not fully account 
for the high-energy behaviour of QCD. 

4. On the basis of Reggeon dominance, we have de¬ 
rived an explicit representation of the subtraction 
function in terms of the electron scattering cross 
sections. The representation requires the asymptotic 
behaviour of the longitudinal cross section to be 
known up to contributions that disappear at high 
energies. For the proton Compton amplitude, where 
Pomeron exchange generates the dominating contri¬ 
bution, the available information does not suffice to 
reliably evaluate the subtraction function. In the dif¬ 
ference between proton and neutron, however, the 
Pomeron drops out. We have shown that the exper¬ 
imental information available at low photon virtual¬ 
ity does suffice to work out the subtraction function 
relevant for this difference. 

5. In [31], the electron cross sections or, ul and the 
subtraction function S\(q 2 ) are instead treated as 
physically independent quantities. The authors in¬ 
voke the low-energy theorem that relates the value 
of the subtraction function at q 2 = 0 to the mag¬ 
netic polarizability and use experimental informa¬ 
tion about the latter to pin down the value of the 
subtraction function at the origin. As direct experi¬ 
mental information about the q 2 -dependence is not 
available, the authors construct a model for that. 
Figure [7] compares their model with our prediction. 
As pointed out in [3B], the model of m is not con¬ 
sistent with the fact that the coefficient of the loga¬ 


rithmic divergence vanishes in the chiral limit. The 
alternative ansatz for the subtraction function pro¬ 
posed there, which does obey this constraint, is also 
shown in figure [7] 

6. The authors of [3TJ use their ansatz for the sub¬ 
traction function to evaluate the difference be¬ 
tween the self-energies of proton and neutron and 
obtain rn^ GM = 1.30(03)(47) MeV, significantly 
higher than the result obtained in [J, m GL = 
0.76(30) MeV. The difference is blamed on a ’techni¬ 
cal oversight’ committed in [5j. This claim is wrong: 
it suffices to replace their ansatz for the subtrac¬ 
tion function with the parameter-free representa¬ 
tion used in [4], which is spelt out explicitly in (l26l) 
above. Leaving all other elements of their calcu¬ 
lation as they are, the central value for the self¬ 
energy difference then drops to m 7 = 0.63 MeV, 
thereby neatly confirming the old result. The same 
conclusion is reached with the calculation performed 
in [TT77I . 

7. We emphasize that the present work only concerns 
low photon virtualities. An update of the analysis 
carried out in [4] which accounts for the progress 
made on the experimental and theoretical sides dur¬ 
ing the last 40 years - in particular an evaluation 
of the contributions from the deep inelastic region 
which accounts for the violations of Bjorken scaling 
- is still missing. 

8. Our representation for the subtraction function also 
leads to a prediction for the difference between the 
electric polarizabilities of proton and neutron. The 
result is given in (l49l) . Using the currently accepted 
results obtained from the Baldin sum rule, this also 
determines the difference of the magnetic polariz¬ 
abilities and, using the comparatively rather precise, 
known value of the electric polarizability of the pro¬ 
ton, we obtain an estimate also for the polarizabili¬ 
ties of the neutron. The result is given in (1571) . 

9. The fact that the results obtained from Reggeon 
dominance are consistent with experiment and even 
somewhat more precise amounts to a nontrivial test 
of the hypothesis that the Compton amplitude is 
free of fixed poles. Quite apart from the possibility 
of taking new data at small photon virtuality, an 
improved representation of the available experimen¬ 
tal information on the cross sections would allow 
us to reduce the uncertainties quite substantially - 
in particular, if the deficiencies of the available pa- 
rameterizations mentioned in subsection 15.11 could 
be removed, the main source of uncertainties in our 
calculation would immediately disappear. 

10. The main problem we are facing with our analysis 
is that all of the well-established features of elec- 
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tron scattering drop out when taking the difference 
between proton and neutron: the leading terms of 
the chiral perturbation series are the same, the con¬ 
tribution from the most prominent resonance, the 
A(1232), is the same, and the leading asymptotic 
term due to Pomeron exchange is also the same. 
Since all of these contributions cancel out, not much 
is left over. Only a fixed pole could prevent the sub¬ 
traction function relevant for the difference between 
proton and neutron from being small. The available 
data do not exclude the occurrence of a fixed pole, 
but they indicate that if the phenomenon occurs at 
all, then the pole must have a rather small residue. 


Since Lorentz invariance, current conservation, and par¬ 
ity only allow two independent tensors of this type, 

, 

jfr = i{( P Y+pv)p-? 

-g^ip-qf-p^q 2 }, (A.4) 

these matrix elements contain two invariants each, 
which only depend on the two variables v = p ■ q/m 
and q 2 (m is the mass of the nucleon). We denote the 
invariants by V\{v, q 2 ), V 2 (v, q 2 ) and Ti(zz, q 2 ), T 2 (i/, q 2 ), 
respectively: 


Acknowledgements We would like to thank S. Alekhin, 
J. Bliimlein, P. Bosted, L. S. Brown, I. Caprini, E. Christy, 
M. Doring, H. Griefihammer, D. Haidt, M. Hilt, G. In- 
gelman, B. Kubis, F.J. Llanes-Estrada, J. McGovern, U.- 
G. Meifiner, C. Merino, O. Nachtmann, D. R. Phillips, 
D. Ronchen, S. Scherer, A. Thomas, L. Tiator, and A. Tk- 
abladze for helpful discussions and correspondence. In par¬ 
ticular, we thank P. Bosted for making the pertinent Fortran 
code available to us. J. G. thanks the Helmholtz-Institut fur 
Strahlen-und Kernphysik at Rheinische Friedrich-Wilhelms- 
Universitat Bonn for hospitality — part of this work was 
performed there. Financial support by BMBF ARCHES, 
the Helmholtz Alliance HA216/EMMI, the Swiss National 
Science Foundation, the Deutsche Forschungsgemeinschaft 
(CRC 16, “Subnuclear Structure of Matter” and CRC 110, 
“Symmetries and the Emergence of Structure in QCD”), 
the Volkswagenstiftung under contract no. 86260, and the 
DOE (Grant No. DE-FG02-00ER41132) is gratefully ac¬ 
knowledged. 


A Notation 

The structure functions are related to the Fourier trans¬ 
form of the spin-averaged matrix element of the current 
commutator, 


V^(p,q) = 2-j d 4 xe^(p\[f(x),fm\p) • (A.l) 

The states are normalized with 


{p',s'\p,s) = 2 p 0 (2tt) 3 6 3 (p' - p )6 s >s (A.2) 

and (p\0\p) stands for y J2 S (P’ s | 0 |p, s )- For the 
Fourier transform of the time-ordered matrix element, 
( p\Tj^(x)j v {y)\p ), we use the normalization 


T^(p,q) = \j 


(A.3) 


V» v (p, q) = Fi(zz, q 2 )K? v + V 2 (v, q 2 )K , 

T^ip, q) = T\ (zz, q 2 )Kr + T 2 (u, q 2 )K^ . (A.5) 


In contrast to the standard structure functions F\, F 2 
the invariants V \, V 2 are free of kinematic singularities 
or zeros. The two sets are related by 


V! = 


—2xFi + F 2 
Fmx 2 v 


V 2 = 


f 2 

2xv 2 


(A. 6 ) 


with x = Q 2 /2mv. 

The notation for the longitudinal structure function 
Fy j is not universal. The convention used in the mini¬ 
review on the structure functions in The Review of Par¬ 
ticle Physics ESI reads 


F L =F 2 - 2 xF x . (A.7) 

The structure functions Viiy.q 2 ) and V 2 {v,q 2 ) rep¬ 
resent linear combinations of the transverse and longi¬ 
tudinal cross sections ctt and ctl: 


Vi = Ni(— Q 2 ctt + v 2 ctl)IQ 2 , V 2 
_ 1 2 mu — Q 2 

1 87 r 2 a em v 2 + Q 2 


N\(cr T + a L ), 
(A. 8 ) 


The value of the structure function V 2 (v, 0) also de¬ 
termines the total cross section for photoproduction, 
0tot = limg2_>.o or: 

CTtot(^) = 47r 2 a em — V 2 (v, 0). (A.9) 


For the one-particle matrix elements of the current, 
we use the notation 


{Pi, si\f\p 2 , s 2 ) = u(pi,si)r M (g)'u(p 2 ,S 2 ), 

r^q) = F D (q 2 )r + F P (q 2 )ia ^, 

2m 

(A.10) 

where q = pi — p 2 . The nucleon spinors are normalized 
with u(p, s')u(p 7 s) = 2m5 s > s . The functions Fp>(t) and 
Fp(t ) are referred to as Dirac and Pauli form factors, 


d 4 xe lq ' x {p\Tf{x)j v ( 0 )|p). 
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respectively. Whenever convenient, we replace these by 
the Sachs form factors, which are defined by 

Gsit) = F D (t) + -^^Fp{t ), G M (t) = F D (t)+F P (t). 

(A.ll) 

In dipole approximation, the form factors are parame¬ 
terized with 

G p E (t) = G d (t), G p M {t) = {l + tf’)G d {t), 

G n E {t) = ^K n G d {t), G n M (t) = n n G d {t ), 

G d (t) = -o , M? ~ 0.71 GeV 2 , (A.12) 

(1 -t/Mj) 2 ' 

where n = Fp( 0) stands for the anomalous magnetic 
moment. 

B Compton scattering 

Virtual Compton scattering in the non-forward direc¬ 
tion provides the bridge between the two processes 
discussed in subsection 14.11 scattering of real photons 
at nonzero scattering angle and scattering of virtual 
photons in the forward direction. Compton scatter¬ 
ing has been thoroughly explored in the literature, for 
the case where both of the two photons are on the 
mass shell (real Compton scattering, RCS) as well as 
when one of them (VCS) or both (VVCS) are off- 
shell [TKlfTiflPITIll^lMlIMlll^) . We normalize the am¬ 
plitude with 

f» v (jpf,Sf,qf\puSi,qi) = l - j d A xe lq f' x 

x {Pf,Sf\Tj^(x)j v (0)\pi,Si), 

(B.13) 

and use matrix notation, collecting the different spin 
orientations in the 2x2 matrix T^(pf,qf\pi 1 qi). The 
spin average is given by the trace of this matrix, 

T^{pf,qf\Pi,qi) = ^ti {T^ipf^flpi^i)} . (B.14) 

B.l Lorentz invariance 

The spin average is not independent of the Lorentz 
frame used. To see why this is so, consider a Lorentz 
transformation: 

U{A)f{x)U~ 1 (A) = (A~ 1 )^ a j a (Ax ). (B.15) 

We denote the pure Lorentz transformation (boost) 
that takes a particle at rest into one of four-momentum 


p by B p and work in the basis where the state \p, s) is 
obtained from the corresponding state at rest, | p, s), by 
application of the relevant boost: \p, s) = U(B p )\p,s). 
Lorentz transformations not only change the momen¬ 
tum, but also subject the spin direction to a rotation, 
referred to as Wigner rotation: 

U(A)\p , s) = Y^ \M s')W s , s {A,p) . (B.16) 

s' 

The Wigner rotation arises because the boost Ba p 
differs from AB p by a rotation, which we denote by 
W(A,p): 

AB p = B Ap W{A,p). (B.17) 

The matrix W s > s (A,p ) in (IB. 161) is the spin tj repre¬ 
sentation of W(A,p). If A is a pure rotation, we have 
W(R,p) = R. Also, since the product of two boosts in 
the same direction is again a boost in that direction, a 
pure Lorentz transformation in the direction of p does 
not generate a Wigner rotation. 

Lorentz invariance implies the transformation law 

A^T^{pf iq/ \p liqi ) 

= W^AprfT^iApf^qf^AqJW^pi). (B.18) 

In the trace, the Wigner rotations only drop out for 
those Lorentz transformations for which W(A,pi) = 
W(A,pf). In general, this condition is violated. Hence 
knowledge of the spin average in one particular frame 
of reference does not in general suffice to determine the 
spin average in a different frame: the transformation 
law (IB. 181) involves the entire matrix T’ A “ / , including 
the spin-flip components of the amplitude. 

In the Breit frame, Pf+Pi =0, the momenta of the 
initial and final states point in opposite directions. The 
boost which takes the Breit frame into the Lab frame, 
where p t = 0, is a pure Lorentz transformation in that 
direction. Hence the change of frame does not gener¬ 
ate a Wigner rotation. Accordingly, the spin average in 
the Lab is determined by the spin average in the Breit 
frame. 

For the Lorentz transformation that takes the Breit 
frame into the centre-of-mass system, however, this is 
not the case: the Wigner rotation generated by this 
transformation for the initial state differs from the one 
relevant for the final state, W(A 1 p i ) ^ W(A,pj). Hence 
knowledge of the spin average in the Breit frame does 
not suffice to evaluate the spin average in the centre-of- 
mass system or vice versa. 
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B.2 Crossing symmetry, parity, and time reversal 

The symmetry of the time-ordered product, 
Tj li (x)j l '(y) = Tj l, (y)j fl (x), implies invariance 
under crossing of the photons! 10 ! 

T^{pf,qf\ Pi ,qi) = T^(p f ,- qi \ Pi ,-q f ). (B.19) 

Invariance under space reflections amounts to: 

T^(pf,qf\pi,qi) = T a ' 3 (npf ,nqf\irpi,Trqi) , 

(B.20) 

where 7 r = diag(l, — 1, — 1, — 1) inverts the sign of 
the space components but leaves the time compo¬ 
nents alone. Time reversal not only inverts the momen¬ 
tum and spin directions, but in addition interchanges 
the initial and final states. Moreover, the amplitudes 
are mapped into their complex conjugate. Exploiting 
the fact that the Hermitian conjugate of the operator 
is given by j v (y)j IJ ‘{x), time reversal invari¬ 
ance can be brought to the form 

T^(.Pf,qf\Pi,qi) 

= 7r^7r^eT a/3 (7rp i ,7T5 i |7rp/,7rgf/) T e _1 . (B.21) 

The superscript T indicates that the transposed ma¬ 
trix is relevant. The matrix e = io 2 flips the spin in 
the initial and final states. For the Pauli matrices, we 
have eer T e _1 = — <x. The above relations lead to the 
following symmetry property of T 1 ' 1 ': 

T^ipf^qflPi^qi) = eT Aiy (p i ,- 9/ |p / ,-g i ) T e~ 1 . 

(B.22) 

It implies that, in the decomposition T^ v = T^l + 
T^ v cr l of the amplitude in the basis spanned by 
1, oq, er 2 , cr 3 , the spin-independent part, T^ v , is even 
under the operation pi -n- pf , q^ —> — qi , qf —> —qf, 
while the spin-dependent part, T)‘ v , is odd. The map¬ 
ping interchanges the Mandelstam variables s = (pi + 
qi) 2 and u = (pi — qf) 2 , but t = (pf — Pi) 2 as well as 
the photon virtualities q 2 , q 2 stay put. While the Breit 
frame is invariant under this operation, neither the Lab 
frame nor the centre-of-mass system have that property. 

B.3 Low-energy expansion 

The elastic intermediate states generate poles in T tlL '. 
In the Mandelstam variables the poles are located at 

10 If the time-ordered product of the currents is replaced by 
the retarded current commutator, crossing symmetry instead 
relates the amplitude to its complex conjugate m¬ 


s = to 2 and u = m 2 . We refer to these contributions as 
Born terms and denote them by Tg', 

T^ = Tg + . (B.23) 

The decomposition is not unique f5%1140111(J2] . The es¬ 
sential property of the Born terms is that they account 
for the elastic singularities. This ensures that the re¬ 
mainder, T , is regular at qt = qf = 0 and can thus 
be expanded in a Taylor series in powers of the photon 
momenta and energies, which the low-energy expansion 
treats as small. The construction described below leads 
to Born terms that are conserved, so that this also holds 
for the remainder l 11 * * ! 

q^T^ = q^T^ = 0. (B.24) 

Keeping P = \{Pi+Pf) fixed, momentum conserva¬ 
tion determines the initial and final nucleon momenta 
in terms of the photon variables However, un¬ 

less P vanishes, energy conservation leads to a nonlin¬ 
ear constraint on the photon energies, so that it is not 
consistent to treat all of these as quantities of O(q). 
The problem disappears for P = 0, i.e. in the Breit 
frame: energy conservation then implies Wj = w/ = to, 
so that the kinematics is unambiguously determined by 
the independent variables which all count as 

quantities of O(q). This is why the Breit frame is the 
preferred frame of reference for the low-energy expan¬ 
sion. For real Compton scattering, a transparent dis¬ 
cussion of this issue is given in 11131 . 

In the Breit frame, the Taylor series in powers of 
the variables uj, q, ; , q^ takes the form: 

= T% V + tT + tT + ... with T^ = 0{q n ). 

(B.25) 

The symmetry (IB.221) implies that the even terms of the 
series Tj, , T 2 ,... are proportional to the unit matrix 
in spin space, while the odd terms T^ v , T 1 -) ,... exclu¬ 
sively contain spin-dependent terms and do not con¬ 
tribute to the spin average. 

11 In general, it is not a trivial matter to impose current con¬ 
servation on the Born terms. If one for instance evaluates the 
one-particle singularities in the space components Tg and 

determines the remaining components of by solving the 
constraints qj) Tg' = q\, Tg' = 0, one in general arrives at 

a representation for the Born terms that contains kinematic 

singularities. The presence of kinematic singularities in the 
Born terms complicates the analysis because the remainder 
T^ is then not regular at qi = qf = 0. Compton scattering 
on the pion illustrates the problem: the one-particle singular¬ 
ities do not generate a term proportional to g , but unless a 
regular term of this type is allowed for, the representation of 
the Born terms can be consistent with current conservation 
only if it contains kinematic singularities. 
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We now turn to the consequences of Lorentz invari¬ 
ance for the low-energy expansion. If the photon ener¬ 
gies and momenta are small of O(q), the momentum 
transferred to the nucleon is small as well. We thus 
only need to consider Lorentz frames where the nucleon 
momenta are also small of O(q). The standard choices 
(Laboratory, centre-of-mass system, Breit frame) all be¬ 
long to this category. If a Lorentz transformation A is 
to connect two such frames, then the relative velocity 
must be small, so that the standard decomposition into 
a boost and a rotation, A = BR only involves a small 
boost. Hence it suffices to analyze the transformation 
properties under rotations and under small boosts. 

The behaviour under rotations is trivial, because the 
corresponding Wigner rotations in the initial and final 
states are identical, W(R,pi) = W(R,pf) = R. Hence 
they leave the spin-independent part of the amplitude 
alone and transform the matrices er = {ay, < 72 , 03 } oc¬ 
curring in the spin-dependent part like a vector. 

To analyze the properties of small boosts, the 
SL(2,C) representation of the Lorentz group is more 
convenient than the one acting on the coordinates and 
momenta. Consider the matrix B = exp(^u • er) g 
SL(2,C), which represents a pure Lorentz transforma¬ 
tion with a small velocity v = O(q). The product of 
two such boosts is given by BB' = l + \w ■ cr + ^w 2 + 
j(v x v') • er + 0(q 3 ), with w = v + v'. At 0(q ), this 
is a pure Lorentz transformation^ but at 0 (q 2 ), the 
product in addition contains a small Wigner rotation: 
W = 1 + |(v X v') • er + 0(q 4 ). 

The essential point here is that, for the pure Lorentz 
transformation needed to remove the relative velocity of 
the two systems, the Wigner rotations generated in the 
initial and final states are at most of order q 2 . Accord¬ 
ingly, the Lorentz invariance condition (IB.181) implies 
that, if A is a boost that only generates velocities of 
O(q), the individual terms of the low-energy expansion 
transform like ordinary tensors - up to higher order 
corrections: 


BaB'pT^f (p f , q f \Pi, q%) = (Bp f , Bq f \B Pi , B Qi ) 

+ 0(q n+2 ). (B.26) 


The corrections only matter if the expansion is taken 
beyond next-to-leading order. In particular, the lead¬ 
ing even term of the series is proportional to the unit 
matrix in spin space even if the reference system is not 
identified with the Breit frame. 


12 The fact that the velocity addition is modified only shows 
up at 0(q 3 ). 


B.4 Spin average 


The Breit frame also offers a convenient decomposition 
of the amplitude into independent tensors. In the fol¬ 
lowing, we explicitly construct the decomposition for 
the spin average, which we denote by . Current 
conservation ( qjj^ T = ql = 0) implies that the 
amplitude is uniquely determined by its space com¬ 
ponents. Furthermore, rotation invariance ensures that 
the spin average involves five independent amplitudes 
(a, b =1,2,3): 

T ab = S ab h+q?q b I 2 +q?q b I 3 +q a f q b h+q a f q b I 5 . (B.27) 

Since this decomposition exclusively makes use of rota¬ 
tion invariance, the coefficients I n are free of kinematic 
singularities. They depend on the rotation-invariant 
quantities u>, |<?j|, |<7/|, and ■ q r , which can be ex¬ 
pressed in terms of the Mandelstam variables and the 
photon virtualities. In view of s + t + u = 2 m 2 + q 2 + 
q 2 , only four of these are independent, for instance: 
s,u,q 2 ,q 2 . 

The trace of the relation (1B.22I) implies that the 
invariants are symmetric under the interchange of s and 

u, 

I n (s, u, q 2 , q 2 ) = I n (u, s, q 2 , q 2 ), n= 1,...,5, (B.28) 

and the crossing symmetry relation (IB. 191) then shows 
that they are symmetric under q 2 <r ■> q 2 as well, except 
that I 3 and I 4 are interchanged, 

In(s,u,<g,q}) = I n (s,u,q 2 f ,q 2 ) , n = 1,2,5, 
h(s,u,q 2 ,q 2 f ) = I 4 {s,u, q 2 f , q 2 ). (B.29) 

The last relation implies that if the photon virtualities 
are the same - in particular for real Compton scattering 
- there are only four independent amplitudes (l3j . 

We add a remark concerning the spin average for 
the case where the scattering amplitude is written in 
the form 


T^{pf,Sf,q f \pi,Si,qi) 

= M(p/,s/)M' i!/ (p / ,g / |p i ,g i )M(p l ,s i ). (B.30) 

If the momenta and pj are parallel, the sum over the 
spins can be represented as a product of projectors: 


^2 u{pi,s) <g) u(p f , s ) = j~, iti + TO ) • (# f + m ). 

S 

K = 2m\]\ — t/Am 2 , (B.31) 


so that spin average can be represented as 

T^ v (pf,qf\Pi,qi) 

= ^ tr {(p f +m) ■ M^(p f ,q f \p i: qi) ■ (fb +m)j . 

(B.32) 
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This formula might suggest that a frame-independent 
definition of the spin average does exist. As already 
noted by Tarrach [15], this is not the case, however: the 
relation (IB.3211 is correct only in those Lorentz frames 
where p i and Pf are parallel. In particular, it does not 
hold in the centre-of-mass frame, where the right hand 
side of (IB.321) does not represent the spin average per¬ 
formed on the left hand side. 


B.5 Born terms 

The Born terms contain poles along the lines s = m 2 
and u = to 2 . The residue of the pole in the s-channel 
involves a sum over the one-particle matrix elements of 
the current: 


U (P/’ S f) r ' 1 (- ( lf) U (.Pn, s)u{p n , s)r U (q i )u{p ll Si) . 

S 

(B.33) 

At the poles, the momentum of the intermediate state 
is on the mass-shell, p 2 = to 2 , and the sum over the 
spin directions is given by 


'^2u(p n ,s)u(p n ,s) = f n +m. (B.34) 


explicit representation for the Born terms: 

R \/l — t/Am 2 
t b — v _ L _ u n — ^ 

(ra z — s){m z — u) 

B x = 4 m 2 u 2 F l D F f D + + F^F l P + F l P F f p ), 

B 2 = w 2 (F i D F f P + F P F P + F P F P ), 

B 3 = u 2 f£(F* d + F ' P ), 

B 4 = w 2 i£(i^ + F f p ), 

^5 = 0 , 

wi = ^(q 2 + q}) 2 - ^t 2 -tw 2 , 

w 2 = \{q 2 i+q))-\t-u 2 . (B.36) 

The symbol F l D stands for F p (q 2 ) and F P ,F P ,F P are 
defined analogously. In the Breit frame, the variable 
ui represents the energy of the photons. Expressed in 
terms of the Mandelstam variables, we have 


s — u 

2y/4 to 2 — t 


(B.37) 


These expressions of course satisfy the relations (IB. 281) 
and (IB.291) . In the case of equal photon virtualities, 
Born terms for the spin-averaged amplitude are also 
provided in [3D]. The representation specified in (3) and 
(4) of that work differs from ours in the overall nor¬ 
malization of the amplitude: basically, it amounts to 
replacing the factor K in (IB.31|) by 2 to. 


We specify the residues off the mass-shell with analytic 
continuation, simply replacing p n with the total mo¬ 
mentum, P = Pi + qi = Pf + qf, also for P 2 ^ to 2 : 


TB'(pf,s f ,qf\pi,s i ,q i ) = 


B^(Pf,s fl q f \pi,q l ) 

TO 2 — S 

B vlx (Pf,Sf, - qi\pi, Si, -q f ) 


m* — u 


1 


B^(pf,s f ,q f \p i ,s l ,q l ) = -u{p f ,Sf)r^{-q f ) 

x (f + m)r v (qi)u(pi, Si). 


(B.35) 


The two terms represent the two tree graphs obtained 
with the standard Feynman rules, except that the 
photon-nucleon vertices are equipped with form factors 
according to (I A. 101) . Taken separately, the two terms do 
not obey current conservation, but taken together they 
do. 

In the Breit frame, the spin average can be evaluated 
with (IB.311) . Comparing the resulting expression for the 
space components with (IB.271) . we obtain the following 


B.6 Leading low-energy constants 


Since u>, q,, q f are independent variables and T 3 U is in¬ 
dependent thereof, the constraint ql f 1 ^ = 0 imme¬ 
diately implies Tg " = 0. The term T x v is a linear 
combination of the variables u>, q i , qj , but this property 
is inconsistent with current conservation: since invari¬ 
ance under space reflections requires the components 
T° x and T ab to be even under a reversal of the pho¬ 
ton momenta, they must be independent thereof. The 
constraint w 2 T°° = q^q b T x b can then only be obeyed if 

T°° and T°^ both vanish, but this is compatible with 
current conservation only if the remaining components, 
T \ a , T x °, also vanish. 

Accordingly, the low-energy expansion of T^ only 
starts at 0(q 2 ). Indeed, the following calculation shows 
that there are exactly two independent conserved ten¬ 
sors of that order. As noted above, the contribution of 
0(q 2 ) is spin independent, T 1 .^ = T$ v 1 and rotation 
invariance requires the space components T% b to be of 
the form (IB.271) with I\ = Ci w 2 +C 2 qf-qi+c 3 q 2 +c 4 q 2 , 
while the other coefficients are constants. For oj = 0, 
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current conservation requires q a fT% b = q\T% b = 0. This 
implies I 2 = — C 2 , I 3 = I 4 = h = C 3 = C4 = 0. For 
the space components, the general solution of the con¬ 
ditions imposed by rotation invariance and current con¬ 
servation at 0 (q 2 ) thus reads 

T“ b = Cl S ab co 2 + c 2 {S ab q f • q t - qU)} ■ (B.38) 

Current conservation then fixes the remaining compo¬ 
nents in terms of the same two constants: 

Tf = Cl Ujq“, f° a = Cl Luq a fl f ° 0 = Cl q r q t . 

(B.39) 

As observed already by Klein [1141 , there are only two 
conserved tensors of polynomial form at 0 (q 2 ). 

B.7 Real Compton scattering 

For real photons, q 2 = q 2 = 0, the projection onto 
the polarization vectors (e° = 0, e ■ q = 0) annihi¬ 
lates the time components as well as the contributions 
from I 3 , 14 , 15 . At low frequencies, the scattering am¬ 
plitude is fully determined by the charge of the par¬ 
ticle m, which we express in units of the proton 
charge, Q e = -Fd(O). At first order in the expansion 
in powers of the photon frequency ui, the anomalous 
magnetic moment k = Fp( 0) also shows up [11611117] , 
and at 0 (ur ), further contributions, characterized by 
the two low-energy constants in (IB.3811 , manifest them¬ 
selves [114 : . In the Breit frame, the low-energy expan¬ 
sion of the spin-averaged amplitude starts with 

= Ci e* ■ Ei + C 2 (e* f ■ m) (e* • n f ), 

Ci = -Ql 

+ ^2 (1 _ -){(1 + z)(Qe + «) 2 - Ql} 

+ W 2 (ci + C 2 z) + O(w^) , 

C 2 = -— yK(2Q e + k)z - oj 2 c 2 + 0(w 4 ), (B.40) 
4 m z 

where n,; = q i /u,nf = qj-/ui are the unit vectors in 
the direction of the initial and final photon momenta, 
respectively, and z = rif ■ rii (in the Breit frame, the 
term linear in ui does not contribute to the spin aver¬ 
age) . Comparison with the well-known low-energy rep¬ 
resentation of the Compton scattering amplitude (see 
for instance (2.1) and (2.5) in [37J) shows that, up to 
normalization, the low-energy constants Ci and c 2 rep¬ 
resent the electric and magnetic polarizabilities, respec¬ 
tively: 

771 m a rr, at \ 

c 1 = - a e, c 2 =- Pm- (B.41) 

^em ^em 


We emphasize that only the sum of the contributions 
from the Born terms and the polarizabilities manifests 
itself in Compton scattering. Both the choice of the ref¬ 
erence frame and the choice of the Born terms are a 
matter of convention [ 38l!40l[T02l . Accordingly, the lit¬ 
erature contains several different variants of the above 
representation. In the context of the present paper, the 
relation (IB.41[) amounts to a definition of the polariz¬ 
abilities. 

Together with the Born terms, the two polarizabil¬ 
ities c/.Ei Pm determine the low-energy expansion not 
only of the spin average but of the entire amplitude 
up to and including 0(q 2 ). In fact, if not only the pho¬ 
ton energies and momenta but the nucleon momenta are 
also booked as small quantities of O(g), this statement 
is valid in any reference frame: the leading terms can 
be written in a manifestly Lorentz invariant manner, 
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=T% + - + ( a E + /3 m)K^} 1 

®em 

+ 0(g 3 ). (B.42) 

The quantities K^ v , K^ represent the generalization of 
the conserved tensors specified in (1A.4I) to nonforward 
directions, 

Kr = QiQf ~ ' Qi > 

K! r = ^{(P IJ Qf + P , 'q l ;)P-q 

- <f v {P ■ qf - P IJ P v q.f -qi], (B.43) 

with P = \(j>i +pf). Conservation of energy and mo¬ 
mentum implies that P ■ qf coincides with P ■ qt = P ■ q. 
While the Breit frame formulae (IB.381) , (IB.391) only con¬ 
tain terms of 0(g 2 ), the representation (IB.421) of the 
contributions from the polarizabilities includes higher 
orders of the low-energy expansion. 


B.8 Low-energy theorems for Tj(i/, q 2 ) and T 2 (u, q 2 ) 

In the forward direction, the spin average is a Lorentz 
invariant notion and only two of the four invariant vari¬ 
ables are independent: the photon virtualities are the 
same, q 2 = q 2 = q 2 1 and the momentum transfer t van¬ 
ishes. The standard variable v = (s — u)/Am coincides 
with the frequency in the Breit frame, ui = v. 

The spin average involves the two invariants Tj, T 2 
defined in (IA.3I) and (IA.5I) . In the notation of (IB.27I) 
these amplitudes are given by 

T\ = I 2 + J3 + I4 + J5 , q 2 T\ + u 2 T 2 = Ii . 


(B.44) 
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While for real Compton scattering only the invariants 
/i, I 2 count, the amplitudes relevant for Compton scat¬ 
tering of virtual photons in the forward direction also 
pick up a contribution from J3, J4,J5. 

The comparison of the above expressions for the 
Born terms with the elastic part of the forward am¬ 
plitudes in (fT5l) shows that the singularities on the left 
and right hand sides of (IB. 441) are indeed the same, but 
in the case of Tf, the regular parts differ: 

Tf = If + if + if + If + ^Fpict 2 ), 

T 2 = ~ 9 2 (4 B + If + If + If)} ■ (B.45) 


The difference also shows up when evaluating the rela¬ 
tion (IB.381) in the forward direction, where it implies a 
low-energy theorem for T\ [1011551I551ll02lll07j as well 
as one for T 2 m, 


Ti nel (0,0) 
Tf e \ 0 , 0 ) 


k m 
Am 2 a em M ’ 

Tfl 

- {ole + /3m) • 


(B.46) 


This demonstrates that, although the amplitudes 
Ti(i/,q 2 ), T 2 (v,q 2 ) do not determine the angular dis¬ 
tribution of real Compton scattering, they do encode 
the polarizabilities. The calculation described in ap¬ 
pendix [BTs] removes the apparent contradiction: the am¬ 
plitudes J3 — If — If , I 5 all disappear if the photon 
energy is set equal to zero. We repeat that the decompo¬ 
sition of the amplitude into a contribution generated by 
the elastic singularities and a remainder only becomes 
unique if the asymptotic behaviour is specified. The ex¬ 
tra term in T\ arises because the recipe used above to 
specify the Born terms implies that the amplitudes if 
do not tend to zero when v —> 00. In the above analy¬ 
sis of non-forward Compton scattering, the asymptotic 
behaviour does not play any role - accordingly, the reg¬ 
ular parts of the Born terms used in that analysis are 
without physical significance. In contrast, the decompo¬ 
sition of the forward amplitudes into an elastic and an 
inelastic part set up in section [3] does invoke the asymp¬ 
totic behaviour. It implies that the polarizabilities do 
pick up a contribution from the elastic singularities. 


C Causality 

The structure functions V\(v,q 2 ),V 2 (v,q 2 ) are experi¬ 
mentally accessible only in the space-like region, q 2 < 0. 
As discussed in detail in [T1 I118I . causality - the fact 
that the current commutator vanishes outside the light- 
cone - very strongly constrains their continuation into 


the time-like region. General properties of causal func¬ 
tions are described in I H). l2()l and explicit represen¬ 
tations that manifestly incorporate causality I121H124] 
are available. These have been used, in particular, in 
the analysis of the structure functions at high ener¬ 
gies [1251 . The further representation constructed in ap- 
pendix[D]shows that the contributions generated by the 
elastic intermediate states or individual resonances can 
be written in manifestly causal form. 

In the present context, the key statement!^ is that 
the continuation is uniquely determined up to a poly¬ 
nomial in the variable v: 

N 

V(v,q 2 ) = ^cr„(g 2 )z/ 2n ’ (C.47) 

71—0 

where the coefficients cr n (q 2 ) vanish for q 2 < 0. In Regge 
language, such contributions represent fixed poles in the 
angular momentum plane, located at integer values of 
the angular momentum. 

Regge asymptotics excludes fixed poles in V 2 (v,q 2 ), 
but a term with n = 0, that is a fixed pole with J = 0, 
is not a priori ruled out in V\(y 1 q 2 ): 

Vl p (v,q 2 ) = e(v)a(q 2 ). (CAS) 


A term of this form is also consistent with the short 
distance properties of QCD, which ensure that, in the 
Bjorken limit, where v and q 2 both become large, the 
structure functions tend to zero - this merely imposes a 
constraint on the asymptotic behaviour, which in par¬ 
ticular requires that cr(s) disappears when s becomes 
large. If a fixed pole were present in Vi(^, g 2 ), it would 
not show up in the electron cross sections, but would 
affect the time-ordered amplitude, through the term 



ds a (s) 
s — q 2 — ie 


(C.49) 


Accordingly, a formula that expresses the electromag¬ 
netic self-energy in terms of the electron cross sections 
could then not be given, nor would it be possible to 
express the polarizabilities of the nucleon in terms of 
these cross sections. 

The analysis of [1] is based on the assumption that 
the matrix element of the current commutator is free 
of fixed poles, so that the electron cross sections unam¬ 
biguously determine the structure functions V\ (v,q 2 ) : 
V 2 (p,q 2 ), not only in the space-like region, but also 
for time-like momenta. The short distance properties 
of QCD ensure that there is then no ambiguity in 
Ti(i/,q 2 ), T 2 (v,q 2 ) either: the electron cross sections 
fully determine these. Accordingly, the electromagnetic 
self-energy as well as the polarizabilities of the nucleon 


13 See theorem 2 in and theorem A in nm 
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are determined by these cross sections, at least in prin¬ 
ciple. 


D A new causal representation 

Consider the product of two retarded propagators 


t re \p,q) = 


Mi - ( [q° + *e ) 2 + q 2 

1 


m! - (P° + 9 ° + *e ) 2 + (p + q ) 2 

The Fourier transform of this amplitude, 


(D.50) 


f et (p,x)= / d\e~™ H re \p,q), 


(D.51) 


is given by the convolution of the two propagators in 
coordinate space. Since these vanish outside the for¬ 
ward light-cone, the same is also true of the convo¬ 
lution: f ret ( p , x ) differs from zero only in the forward 
light-cone, x 2 > 0, a : 0 > 0. 

This reflects the properties of the integrand 
in (ID.bill , which contains four poles 

f et (n a) =_-_ 

[P,q) ((jO-wi+teXgO + Wi+ie) 

1 

(q° +p° - uj 2 + ie) (q° + p° + oj 2 + ie) ’ 

wi = \Jp\ + q 2 , w 2 = yjt4 + (q + p) 2 • 

(D.52) 

All of these occur in the lower half of the ( 7 °-plane. 
The path of integration can therefore be deformed into 
a segment from —oo to —R, a semi-circle of radius R 
and a segment from R to + 00 . If x° < 0, the factor 
e -iq x su pp resses the integrand if R is taken large, 
except for the segments where the imaginary part of q° 
is not large. But there, t ret (p, q) is small, of order 1 /R A . 
Since the length of these segments is of order R, their 
contributions also tend to zero if R is taken large. Since 
the integral is path-independent, it vanishes for x° < 0 . 

The quantity t let (p,x) is Lorentz invariant. The 
time component of the vector x depends on the frame 
chosen, but for any point outside the forward light-cone, 
there is a frame where the time component is negative. 
This confirms that f ret (p, x) vanishes outside the for¬ 
ward light-cone. 

The advanced version of the amplitude only differs 
in the sign of the ie prescription. Like the advanced 
propagators in coordinate space, the Fourier transform 
of f adv (p, q) is different from zero only in the backward 


light-cone, x 2 > 0, x° < 0. The difference between the 
two, 


V(p, q) = f et (p , q) - t adv (p, q)) , (D.53) 

2 m 

is therefore causal: the Fourier transform of v(p , q) van¬ 
ishes outside the light-cone. Since t adv (p, 9 ) is the com¬ 
plex conjugate of t ret (p,q ), the function v(p : q) is real. 

For space-like momenta, v(p, q) picks up a contribu¬ 
tion only from the poles at q° = —p° ± uj 2 : 


v(p,q) 


2 _ 2 S((p + q ) 2 - P 2)1 9 2 <0,i/>0. 

Mi q 

(D.54) 


This demonstrates that any function which for space¬ 
like momenta can be represented as 


/* OO f>C 

V(",q 2 ) = da 
Jo Jo 


db 


P(a,b) 

a — q 2 

x {S(q 2 + 2 vm — b) — 5(q 2 


2 vm — 6 )} , 

(D.55) 


admits a causal continuation into the time-like region 
(we have replaced p\ and by a = p 2 , b = p 2 ~ m2 > 
respectively and imposed the condition b > 0 , which 
ensures that, for space-like momenta, the support of 
V(v,q 2 ) is contained in the physical region, Q 2 < 
2m|^|). There is only one continuation that is free of 
fixed poles. The corresponding time-ordered amplitude 
is given by 


/•OO f- OO 

T(v, q 2 ) = da db 

Jo Jo 


P(a,b) 


0 a — q“ — ie 

1 


b — q 2 + 2 vm — ie 

+ 1 ) 

b — q 2 — 2 vm — ie [ 


(D.56) 


This representation, in particular, yields a causal de¬ 
scription of the elastic contributions discussed in sec¬ 
tion [3] Since the form factors are analytic functions of 
t, the coefficients uf 1 , vf introduced in (fTOl) do admit a 
representation of the form 

vf(q 2 ) = [ da — . , (D.57) 

Jo a- q 2 -ie 

so that the spectral functions pi(a,b) = ai(a)S(b) in¬ 
deed generate the elastic contributions. 
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E Errata in the literature 

In the present appendix, we rectify a number of incor¬ 
rect statements made in the literature [301133] about the 
work of Gasser and Leutwyler 0. Unfortunately, some 
of these already propagated EMM®. 

The only deficiency of the analysis reported in 0 
which we are aware of concerns the evaluation of the 
contributions arising from the deep inelastic region: the 
violations of Bjorken scaling are not accounted for. In 
particular, if the ratio x = Q 2 /2mv is kept fixed, the 
longitudinal structure function Fl = F 2 — 2xF\ is as¬ 
sumed to tend to zero at high energies, in inverse pro¬ 
portion to Q 2 . In QCD, Fl only disappears in propor¬ 
tion to a 3 oc l/log(<5 2 ). The calculation yet needs to 
be improved to account for the scaling violations, but 
we doubt that this will significantly affect the numerics 
in 0 , the contribution from the entire deep inelastic 
region was found to be small. 

The essence of the analysis in 0 is recapitulated in 
section [21 Reggeon dominance implies that the subtrac¬ 
tion function is uniquely determined by the cross sec¬ 
tion of the reaction e+p —> e + X. In [30], for instance, 
this crucial point is not addressed at all. Instead, the 
paper contains the following statement: 

In their work, they acknowledged the need for a sub¬ 
tracted dispersion relation but proceeded to ignore this 
issue, as the subtraction constant could not be com¬ 
puted. 

Quite the contrary, the analysis in 0 not only includes 
an explicit evaluation of the subtraction occurring in 
the fixed-g 2 dispersion relation for Ti, but contains a 
detailed discussion of the matter in a separate section, 
entitled “Regge poles, fixed poles, subtractions and the 
like.” 

In [3T] the authors write: 

In 0 it was claimed that the elastic contributions to 
1 1 could be evaluated with an unsubtracted dispersive 
analysis. However, performing an unsubtracted disper¬ 
sive analysis of the elastic contributions [... ] leads to 
inconsistent results. 

This conclusion is obtained by comparing the decompo¬ 
sition of T^ v fp, q) specified in appendix lAl CTi = —\t\, 
T 2 = \t'i in their notation) with an alternative decom¬ 
position, introduced by the authors as Ti = 2q 2 Ti + 
2v 2 T 2 , T 2 = —2 q 2 T 2 . While the amplitudes T\,T 2 are 
free of kinematic singularities and zeros, Ti, T 2 are not. 
That is why the authors run into an inconsistency: 
while T 2 does obey an unsubtracted dispersion rela¬ 
tion, the dispersive representation for v 2 T 2 requires a 
subtraction. The need for a subtraction also shows up 
in the asymptotic behaviour of the elastic contributions 
in (El), which implies that does not disappear when 


v — > 00 and can therefore not possibly obey an unsub¬ 
tracted dispersion relation. The calculation described 
by the authors merely shows that amplitudes with kine¬ 
matic zeros may fail to obey unsubtracted dispersion 
relations. It does not demonstrate that the represen¬ 
tation of the elastic contributions Tf 1 , Tff in 0 is in¬ 
correct - in fact, as shown in appendix [Dj even taken 
by itself, that representation for the elastic contribu¬ 
tions to Ti(v,q 2 ), T 2 {v,q 2 ) obeys all of the constraints 
imposed by causality. 

Further, in the conclusion of m, the authors state: 
A technical oversight in the evaluation of the elastic 
contribution was highlighted resulting in a larger cen¬ 
tral value than previously obtained 0 . 

The claim of a ’technical oversight’ suggests that 
the value 0.76 MeV for the elastic contribution ob¬ 
tained in 0 is incorrect and needs to be replaced by 
1.39(2) MeV. That, however, is not the case. The num¬ 
ber quoted in 0 concerns the full contribution of the 
elastic intermediate states to the Cottingham formula, 
which includes the elastic contribution to the subtrac¬ 
tion function. In [31], the elastic part of the mass-shift 
is instead identified with the contribution from the elas¬ 
tic states to the subtracted dispersion integral. In fact, 
the authors notice that their model for the subtraction 
function also contains elastic contributions and that if 
these are accounted for, the claimed discrepancy in the 
value of the elastic contribution to the self-energy differ¬ 
ence disappears 03 This demonstrates that their claim 
is wrong and that the evaluation in 0 is correct. 

The incorrect statements concerning the analysis 
in 0 are iterated elsewhere. For instance, the following 
claims are made: 

[32] : One is lead to conclude that the unknown subtrac¬ 
tion function can not be evaded and the evaluation of 
8mA in Ref. 0 is not correct. 

[33] : This work uncovered a technical oversight in 
work of Gasser and Leutwyler 0£D related to a sub¬ 
tracted dispersion integral, which unfortunately invali¬ 
dates their result. 

Finally we point to an erratum in [32] : 

In Ref. 0 an argument to evade the subtraction func¬ 
tion based on the parton model was presented. However, 
as was first noted in Ref. [21], the argument was based 
on false assumptions about the scaling violations of the 
Callan-Gross relation m- 

The necessity of a subtraction arises from the Regge 
behaviour at fixed q 2 . The analysis in 0 relies on stan- 

14 Incidentally, although their amended expression for the 
elastic part of the self-energy is numerically close to the con¬ 
tribution generated by the elastic singularities, this is only 
approximately so: while these singularities are proportional 
to Gg and G\j, their expression for the elastic part of the 
mass-shift in addition involves an integral over GeGm- 
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dard Regge behaviour and hence necessarily involves a 
subtracted dispersion relation for T\. Scaling violations 
concern the behaviour of the structure functions in the 
deep inelastic region, where q 2 becomes large. It is true 
that in [4], the scaling violations are ignored, but to 
claim that this is used in [3] as an argument to evade the 
subtraction function is plain wrong: a subtraction was 
made and the subtraction function was determined. 
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